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ABSTRACT 

We have used a combination of broad-band near-infrared and optical Johnson-Cousins photometry to study 
the dust properties in the line of sight to the Galactic globular cluster M4. We have investigated the reddening 
effects in terms of absolute strength and variation across the cluster field, as well as the shape of the reddening 
law defined by the type of dust. All three aspects had been poorly defined for this system and, consequently, 
there has been controversy about the absolute distance to this globular cluster, which is closest to the Sun. Here, 
we determine the ratio of absolute to selective extinction {Ry) in the line of sight towards M4, which is known to 
be a useful indicator for the type of dust and therefore characterizes the applicable reddening law. Our method 
is independent of age assumptions and appears to be significantly more precise and accurate than previous 
approaches. We obtain A (//^(B - V) - 3.76 ±0.07 (random error) for the dust in the line of sight to M4 for our 



set of filters. That corresponds to a dust-type parameter Ry - 3.62 + 0.07 in the Cardelli et al. ( |1989 i) reddening 
law. With this value, the distance to M4 is found to be L80±0.05 kpc, corresponding to a true distance modulus 
of (m - M)o - 11.28 ± 0.06 (random error). A reddening map for M4 has been created, which reveals a spatial 
differential reddening of SE(B - V) > 0.2 mag across the field within 10' around the cluster centre; this is 
about 50% of the total mean reddening, which we have determined to be E(B - V) = 0.37 + 0.01. In order to 
provide accurate zero points for the extinction coefficients of our photometric filters, we investigated the impact 
of stellar parameters such as temperature, surface gravity and metallicity on the extinction properties and the 
necessary corrections in different bandpasses. Using both synthetic ATLAS9 spectra and observed spectral 
energy distributions, we found similar sized effects for the range of temperature and surface gravity typical 
of globular cluster stars: each causes a change of about 3% in the necessary correction factor for each filter 
combination. Interestingly, variations in the metallicity cause effects of the same order when the assumed value 
is changed from the solar metallicity ([Fe/H] - 0.0) to [Fe/H]=-2.5. The systematic differences between the 
reddening corrections for a typical main-sequence turnoff star in a metal poor globular cluster and a Vega-like 
star are even stronger (~ 5%). We compared the results from synthetic spectra to those obtained with observed 
spectral energy distributions and found significant differences for temperatures lower than ~ 5 000 K. We have 
attributed these discrepancies to the inadequate treatment of some molecular bands in the B filter within the 
ATLAS9 models. Fortunately, these differences do not affect the principal astrophysical conclusions in this 
study, which are based on stars hotter than 5 000 K. From our calculations, we provide extinction zero points 
for Johnson-Cousins and 2MASS filters, spanning a wide range of stellar parameters and dust types. These 
extinction tables are suited for accurate, object-specific extinction corrections. 

Subject headings: globular clusters: individual (M4, NGC 6723) — ISM: dust, extinction — Techniques: 
photometric 



1. INTRODUCTION 

Globular cluster systems rank among the oldest objects 
known in our Galaxy. Their ages set a lower limit for the 
age of the Universe. Unfortunately, a precise measurement of 
the absolute age of a globular cluster (GC) is strongly depen- 
dent on the knowledge o f its precise dista nce and reddening 
( |Gratton et a l. 2003 , Bolt e & Hogan|1995| l. 

M4 is a peculiar cluster in terms of .nterstellar redden- 
ing. It has a su rprisingly la rge total amount of reddening 
(E(B-V) = 0.35, |Harris|l996l l for its small distance (~ 2 kpc), 
which is due to its location in the Galactic plane, behind 
the Sco-Oph cloud complex. Moreover, the cluster suffers 
from a significant amount of spatially differential redden- 



ing (Cudworth & Rees 



19901 iprake et al.|fl994; and [Ivans 
lEr 



|et al.||1999[ ) which has been taEn into account i n some re- 
cent studies of M4 (e.g. ' 



Marino et al..2008,.Mucciarelli et al. 



201 1 1. In the literature, approximations for peak-to-peak dif- 



ferences within a distance to the cluster centre of at least 10' 



range from 6EiB - V) > 0.05 in Cudworth & Rees ( 1990 1 to 



6E(B - y) = 0.25 in| Mucciarelli et al.| ( [20n] l. 

The extinction law that is used to correct for interstellar red- 
dening can be different for each line of sight, depending on 
the type of dust located between the observer and the object, 
where the total to selective extinction Ry = Av/E(B - V) is 
known to be a goo d index of the extin ction properties of a 
paiticulai- dust type ( [Cardelli et al.|1988| jMathis 1990). How- 
ever, in many studies of star clusters a standard reddening law 
is applied, assuming dust characteristics typical of the diffuse 
interstellar medium (ISM) as the main origin of interstellar 
extinction: a reddening law with a value of Ry close to 3.1 
is typically adopted (e.g. Sneden et al. 1978, McCall 2004). 
However, since Ry defines both the shape of the reddening 
law and the ratio between extinction and reddening, a recon- 
sideration of this standard assumption is appropriate when de- 
termining the absolute distance of a reddened stellar system. 

For M4, there are several hints in former studies where the 
authors suggest an abnormal dust type, or at least remark on 
certain discrepancies that arise when using a standard assump- 
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tion of Ry = 3.1. For example, Ivans et al. ( 1999)1 compare 
their spectroscopically derived temperatures of bright red gi- 
ant branch (RGB) stars to the temperatures derive d from pho- 
tometric ind ices and found Ry - 3.4 ± 0.4. Dixo n & Long-| 
more| ( 1993 ) propose an Ry "closer to 4 then to 3" by evaluat- 
ing the relative location of RGB sequences fr om M4 com- 
pared to M3 
(|1983||. Interestingly. 



Ml 3 and M92 using data from [Frogel et al. 
the only study that does not rely on 
the ca hbration of photom etric magnitudes reaches a similar 
result: [Peterson et al.| ( [T99 5 ) derived a geometric distance for 
M4 from proper motion and radial velocity measurements and 
find a value significantly smaller than the accepted photomet- 
ric distance, a result which is only in agreement with canoni- 
cal horizontal branch (HB) magnitude measurements if a red- 
dening law with Rv ~ 4-is assumed. In addition to direct mea- 
surements of Ry using cluster members, there are several indi- 
rect measurements using individual stars in the Sco-Oph dust 
cloud complex with a line of sigh t close to M4, which yield 
values of Ry aro und 4 (see e.g 
|Vrbaetal.|1993] l. 

Indeed, some studies of M4's colour-magnitude diagram 
(CMP) already use a non-standard, higher value of Ry. 



Clayton & Cardelli 1988 



Richer et al. ( 1997 ) (and later Richer et al. 2004 ) were among 
the first to use a value of Ry - 3.8 for their distance estimate 
when determining the age of M4 from the white dwarf cooling 
sequence. Ho wever, this choice is s olely based on the vague 
statements of 'Peters on et al.| ( |T995] l, and |Vrba et al.|(|1993|l 
who only estimate the value to be "ar ound 4." Later, |Hansen| 
(et al.| (2004 ) and Bedin et alT] ( |2009| ) follow the argumenta- 
tion of Richer and additionally use the results from C layton] 
|& Card elli ( 1988 ), who found a value of Ry = 3.8 for a star 
only one degree away from the li ne of sight of M4 to jus- 
tify their choice. However, Mathis fT990 1 warns in his review 
about the properties of interstellar dust that "it is not possible 
to estimate Ry quantitatively from the environment of a line 
of sight," since significant variations in the properties of dust 



can occur even on these small angular scales (see also Vrba & 
Rydgren|1985 i. Moreover, this value for Ry is based on the 



measurement of only one star, using obsolete stellar models 
to estimate its atmospheric parameters. 

M4 is the closest GC to the sun. Its closeness and low de- 
gree of crowding make the stellar population accessible to es- 
pecially deep photometric and spectroscopic analyses. For 
example, M4 hosts one of the largest populations of identi- 
fied white dwarfs (WDs), making it attractive for absolute age 
determinations and an excellent la boratory f or testing stellar 
evolutionary theory (e.g. Richer et al.[199 7)). It has a signifi- 
cant population of RR Lyrae stars and it shows a bimodal HB, 
with well populated blue and red parts, making it an interest- 
ing object to study the so called "Sec ond Parameter Problem" 
(see e.g. the review by Catelan|2009| or some recent work by 
[Marino et al.|20n] ). 

There is clearly a need for a detailed investigation of the ab- 
solute reddening and the specific reddening law for M4, pro- 
viding more precise information about its dust type and the 
effects of spatial differential reddening across the cluster face. 
A better determination of the absolute photometric properties 
of M4's stars — and consequently its distance — will increase 
the cluster's status as a most attractive target for astrophysical 
investigations. 

In this study we investigate the type of dust in the line of 
sight to M4, as characterized by the appropriate value of Ry 
assuming the validity of the general reddening law given in 
CardelU et al.[ (| 1989| l. We achieve high quantitative preci- 



sion by using a combination of near infrared (NIR) J and Ks 
and optical Johnson-Cousins UBV{Rr)c photometry for the 
cluster M4, supplemented with equivalent photometry for the 
cluster NGC 6723 with which we verify our models in unred- 
dened condition s. We use the n ewest set of Victoria-Regina 
isochrones (VandenBerg et al.j2012 ) to determine colour ex- 
cesses in filter combinations that are sensitive to the type of 
dust and consequently to the shape of the reddening law. To 
increase the precision of parameters such as distance, age and 
the absolute locus of the observed fiducial sequence, we map 
the variations in the reddening across the field of M4 that has 
been surveyed and then correct for those differential effects. 

Since extinction is a function of wavelength, the accurate 
extinction properties of our broadband filters depend upon 
the distribution of observed stellar flux within the passband. 
To minimize systematic errors in our extinction corrections 
and to the determined value of Ry, we investigate differential 
changes in the extinction properties of optical and NIR fil- 
ters and filter combinations appropriate for stars with different 
intrinsic spectral-energy distributions as determined by their 
temperature, surface gravity and metallicity. Specifically, we 
consider whether a single reddening law is sufficient to deter- 
mine cluster parameters from photometry, or whether a star- 
by-star correction is required to avoid significant systematic 
errors. These results are then used to define an object-specific 
reddening law for M4 where the correction zero points are 
tailored for the atmospheric and chemical parameters of this 
cluster instead of assuming a Vega-like star, as is the case for 
most literature values. To test the consistency of the results 
obtained with synthetic and observed stellar fluxes we com- 
pare synthet ic spectra from the ATLAS9 library (Castelli &] 
'Kurucz'2003 ) to observati onal databases from |Pickles , ( ,199^ 
and Sanchez-Blazquez et al. (2006). The outcome of such a 
comparison not only constrains the validity of our results, but 
also reveals possible weak spots in the atmospheric models. 

In ^ we summarize the data and the reduction process to- 
gether with the criteria for choosing our high-quality photo- 
metric sample. Our procedure for generating a reddening map 
for M4 and the spatial differential reddening corrections is ex- 
plained in ^ while in ^ we investigate the effects of stellar 
temperature, surface gravity, metallicity and extinction on the 
shape of the reddening law and calculate an object-specific 
law tailored for M4. Finally, in 0] we characterize the dust 
type toward M4, derive the cluster s absolute distance, and 
discuss random and systematic errors in detail. A summary 
is given in ^ Appendices provide fiducial sequences for M4 
and NGC 6723, and explicit extinction tables for a wide range 
of Teff, [Fe/H] and Ry. 

2. DATA AND DATA REDUCTION 
2.1. Instrumental Photometry 
For this work, we use near infrared J and KV\ ground- 



based photometry for the globular clusters M4 (NGC 6121) 
and NGC 6723. The images for M4 were taken in 2002 with 
the large and small field of SOFI, the infrared imaging camera 
on the New Technology Telescope of ESO on Cerro La Silla 
in Chile as part of the program 69.D-0604(A) (PI: Zoccali) 
and cover a field of ~ 13' x 13'. For NGC 6723 the same 
telescope and instrumentation were used in 2005 as part of 
the program 075.D-0372 (PI: Ferraro). The pre-processing of 
the CCD images for both clusters included bad pixel masking, 

' Throughout the remainder of this paper K and Kg are u.sed interchange- 
ably where both refer to the 2MASS bandpass. 
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bias and dark frame subtraction, and flat-field corrections, as 
well as sky subtraction. 

We use the data reduction packages DAOPHOT II ("Stet-I 
[son 1987; Stetson & Harris 1988) and ALLFRAME (Stetson 
[1994) and related routines to obtain the instrumental photom- 
etry for eac h cluster and transform it to the 2MASS photomet- 
ric system ( [Skrutskie et al.|[2006| l. Since the data reduction 
process for our two clusters is very similar, we will describe 
this work in detail only for M4 as an example. 

In total, 50 frames in each filter were taken with a main 
concentration on the central region of the cluster, which was 
observed 10 times, and one region outside the core, which was 
observed 20 times. The rest of the cluster is covered by only 
one exposure in each filter with slight overlap. 

For the profile-fitting photometry we use a point spread 
function (PSF) which varies quadratically with the position 
in the frame because the SOFI large-field camera was not 
properly aligned, so that the images are radially elongated 
in the left part of each frame. This distortion affects a verti- 
cal strip of about 150 pixels and smoothly disappears toward 
the centre of the frame. Later, we correct the residual spa- 
tial inhomogeneities to first order by applying a differential 
reddening correction (see section [3]). The PSF is calculated 
using a sample of bright, isolatedprogram stars, uniformly 
distributed over the individual frames. Although we have 
written a pipeline to process all frames through DAOPHOT 
and ALLFRAME automatically, we pick out the PSF stars for 
each frame by eye to assure their quality and rule out blend- 
ing. The same stars later serve as local standards to calibrate 
our instrumental photometry to the 2MASS standard system. 

We use DAOMATCH and DAOMASTER to find positional 
transformation equations between the frames, accounting for 
differences in scale and rotation as well as transverse offsets. 
Those equations are then used to match up stars between dif- 
ferent frames in order to get a complete star list for our cata- 
log. 

To calibrate the frames among themselves, DAOGROW 
(Stetson 1990 1 is used to obtain the total integrated instru- 
mental magnitudes (i.e., the total photon count) for our stan- 
dard stars by multiple aperture photometry. These total mag- 
nitudes are then calibrated by direct comparison to the pub- 
lished 2MASS magnitudes for the same stars, as described in 
the next section. 

2.2. Calibration to the 2MASS Standard System 

Our instrumental photometry for M4 must be calibrated 
to the 2MASS standard system: for accurate comparisons 
of our observed data to stellar evolutionary models, it is es- 
pecially important that systematic deviations between instru- 
mental and standard magnitudes are small. 

The same hand-selected, bright and isolated program stars 
used to determine the PSF above now serve as local stan- 
dards for our photometry. As an additional criterion, we 
only use stars below the 1% linearity limit of the detector 
(< 13 000 ADU), and only those that could be cross-identified 
with a 2MASS standard star with a coordinate deviation less 
than 0.3". For M4, these comprise about 360 stars covering 
the whole observation field. Since they are distributed homo- 
geneously over the different frames, we are able to detect and 
compensate for photometric differences between different ex- 
posures as well as spatial inhomogeneities within individual 
fields. 

DAOMASTER is then used to match up 2MASS stars with 
corresponding local standards in our photometry. The equa- 
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Fig. 1. — Photometric differences between our photometry for M4 and 
2MASS standard stars where AJ = J - Jim AS s ™d AK = K- Kimas s ■ Only 
stars are shown that could be cross-identified with a 2MASS star to within 
0.3". The solid black line indicates a theoretical zero offset. Stars used as 
local standards are viewed in black, grey dots indicate all other identified 
2MASS stars from our final high-quality sample (see following paragraphs). 
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Fig. 2. — As in the previous plot, except that photometric differences are 
plotted as a function of V - A" to examine whether there are any trends with 
temperature. 



tion that we use to transform the instrumental magnitudes (m) 
to the standard system (M) is of the general form 



mi = Ml + Zc + acX + Q.(mi — TO2) 



(1) 



and includes a zero-point correction (zc), a term for atmo- 
spheric extinction (a^; with airmass X), and a colour-sensitive 
term (q) to account for differences in filter characteristics be- 
tween the two systems. For our J photometry we find Zc - 
1.956 ± 0.004, = 0.016 ± 0.034, and c, = -0.092 + 0.005. 
For the K filter, Zc = 2.641 + 0.006, a, = 0.126 ± 0.046, 
and Cc = 0.008 + 0.006. The correction coefficients for at- 
mospheric extinction (a^) are a little low for J and quiet high 
for K, but since the variation in airmass within our individ- 
ual frames is only 0.028 for both NIR bandpasses, the im- 
pact on our photometry compared to "standard" assumptions 
of flf a; 0.04 for each filter would only be in the order of 
~ 0.001 mag and therefore negligible for our results. Once 
the transformation equations to the standard system are de- 
fined for each filter with the standard star sample, they are 
applied to all of the stars in the cluster. 

A comparison between our photometry and the 2MASS 
standard system is shown in Figures [T] and [2] where the de- 
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viation between the two systems is plotted as a function of 
magnitude and colour. We find neither a significant zero-point 
offset between the two systems nor a clear trend with colour 
or magnitude, which would arise from insufficient correction 
of the different filter characteristics. 

2.3. Optical Photometry 

The optical Johnson-Cousins photometry for M4 is a com- 
pilation of 4794, primarily archival, images from 48 inde- 
pendent nights distributed among 11 different observing runs. 
The images were obtained between 1994 and 2007 with sev- 
eral different ground-based telescopes ranging from 0.9 m to 
3.6 m, and were analyzed and calibrated as described in Stet^ 



son! (2000 1 and Stetson (2005). In total, the optical database 
covers a held of 25' x l5' around the cluster centre, which 
is about four times the area that we cover with our NIR pho- 
tometry. The database consists of ~ 81 000 stars and the pho- 
tometry reaches down to B ~ 25 mag, V ~ 23.5 mag, and 
/ ~ 21 mag. We use DAOMASTER to match up stars be- 
tween the two photometric sets, where only stars are consid- 
ered for which the coordinate deviation between the optical 
and the NIR database is smaller than 0.3". 

2.4. Selection of a High-Quality Sample 

The initial output of the photometry packages DAOPHOT 
and ALLFRAME yields a total of ~ 21 000 stars with pho- 
tometric information in both J and Ks and in at least two of 
the optical bands. This sample still includes stars with high 
photometric errors, field stars in the line of sight to the actual 
GC system, and variable stars. Since we are more interested 
in photometric precision than in completeness, a high-quality 
subsample (hereafter: HQ-sample) was selected taking into 
consideration the aspects mentioned above. 

To select only stars with the best photometric standard er- 
rors, we have determined the 20-th percentile of standard er- 
rors in magnitude bins of ~ 1 mag straight line segments were 
then interpolated between these points to define a rejection 
criterion as a function of magnitude. This has been done for 
each filter individually to account for the different character- 
istics of each bandpass and the different observing conditions 
for each dataset. It is important to note that we do not auto- 
matically choose a lower boundary for brighter stars. A less 
strict boundary for HB stars has been applied so as to obtain 
a larger, more complete sample in this evolutionary stage and 
thereby to better to define the observed zero-age horizontal 
branch (ZAHB) luminosity, which is an important distance 
and age indicator for the cluster. Furthermore, some pixels in 
the brightest stars occasionally exceed the saturation limit of 
the detector, which increases their observed photometric stan- 
dard error. To produce a clear RGB sequence extending as 
far as possible towards brighter magnitudes, we applied a less 
strict photometric rejection limit for these stars, as well. 

Since the distribution of frames over the cluster field has an 
important impact on the final photometric standard error of a 
star — a large number of frames per observing field decreases 
the photometric error of stars therein — the final HQ-sample 
is not distributed homogeneously over the field, but favours 
areas that have been observed most frequently (see Figure[3]). 

A significant fraction of field stars can be a problem in the 
analysis of GC CMDs, since they are not located on the evolu- 
tionary sequence of the population. Therefore, they can smear 
the actual sequence and, in the worst case, lead to misinterpre- 
tations of observed features. By investigating our CMDs of 
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Fig. 3. — Location of stars from the HQ-sample in the field of M4: The final 
selection (black) is not a homogeneous distribution over the whole cluster 
field (grey), but shows a preference for stars lying in areas with the most 
exposures. 

the whole cluster sample and comparing the total number of 
stars to those which clearly do not follow the common cluster 
sequence, it is clear that the fraction of field stars lying in the 
line of sight to the object is very small for the field of view 
of our observations. From this ratio, we infer the fraction of 
field stars to be < 1%. Therefore we assume that field stars do 
not affect our photometric analysis. 

RR Lyrae stars have been identified by their coor- 
dinates from t he most recen tly updated database of C. 
Clement ( jClement et al. 2001 2011 edition). Only cross- 
identifications with a coordinate deviation of less than 1" and 
only detections lying ~ +1 mag around the ZAHB V-band 
magnitude have been considered to be true RR Lyrae identi- 
fications. By these criteria, 35 RR Lyrae stars were identified 
in our sample, which are plotted in Figure |4] as black circles. 
Except for a few outliers, these stars fall very well in the pre- 
dicted instability strip. For all RR Lyrae stars, we adopt the 
ALLFRAME magnitudes from our photometry and not the 
mean magnitudes. 

In our final HQ-sample, the photometry reaches to well be- 
low the main sequence turn-off (MSTO) to a magnitude of 
~ 22.5 in V and ~ 18.5 in J. Notably, it reaches below the MS 
knee visible only in NIR photometry which is an interesting 
feature for the determinat ion of the distance s, and hence ages, 
of globular clusters (see Bono et al.|2010[ l. At the luminous 
end, the photometry is limited by the linearity limit of the NIR 
data which lies at ~ 11.5 mag in V and ~ 9.0 mag for the J 
filter. 

3. DIFFERENTIAL REDDENING 

In our photometry of M4, both the HB and the RGB se- 
quences are poorly defined (seen as a large scatter in magni- 
tude on the HB and a large scatter in colour along the RCjB) 
compared to the well defined subgiant branch (SGB) and MS. 
Additionally, in CMDs with a combination of optical and NIR 
filters, the colour scatter of the MS increases with luminos- 
ity towards the TO. Considering only the photometric uncer- 
tainty, one would expect the result to be the opposite. Fur- 
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Fig. 4. — Different CMDs containing the HQ-sample of M4. Black circles 
indicate the identified RR Lyrae variable stars in the sample. In the bottom 
plot, the MS "knee," which is seen only in the NIR, is just visible at / ~ 
17.5 mag. The apparent larger spread in colour showed by the NIR filter 
combination is only caused by the smaller colour range covered by this CMD. 
The black bars indicate the length and direction of the reddening vector for 
each filter combination assuming E{B - V) = 0.37 and a standard reddening 
law with i?v = 3.1. 



thermore, the scatter on the HB is much stronger in V than 
it is in the 7-band, suggesting that those effects are mainly 
caused by differential reddening instead of photometric er- 
rors. Generally, differential reddening induces the stars to 
scatter much more around the actual sequence than would be 
expected from their photometric uncertainty. This additional 
scatter, however, is not uniformly evident throughout the GC 
fiducial sequence. In fact, it depends on the angle between 
the sequence and the reddening vector. Regions of the CMD 
where this angle is large show a larger scatter than regions 
where the sequence and the reddening vector are almost par- 
allel. The direction of the reddening vector depends upon the 
actual filters that are used and the extinction relation between 
these filters, defined by the reddening law. In Figure |4] we 
show the reddening vector for each individual filter combina- 
tion using a standard reddening law with TJy = 3.1 for purpose 
of illustration. Note, that for M4 we expect a slightly different 
orientation of the reddening vector. 



3.1. Correction Procedure 

The general idea for determining and correcting spatially 
differential reddening across the face of the cluster is to calcu- 
late for a sample of reference stars the distance from a fiducial 
sequence in a CMD as measured along the reddening vector. 
In theory, this displacement depends only upon the individual 
reddening of the star and is therefore directly correlated with 
its spatial location in the cluster field. Now, the local redden- 
ing value defined by these reference stars can be provisionally 
applied t o all program stars in their immediate ne ighbourhood 
(see e.g., Piotto et al. 1999 Sarajedini et al.|200 7|l. Such a cor- 
rection for spatially ditt'erential efi'ects in the photometry po- 
tentially includes a first-order correction for poor flat-fielding 
or spatially varying PSFs within different frames as well, in 
case they were not modelled perfectly during the data reduc- 
tion process (as is likely the case for our photometry). 

Our strategy for correcting differential reddening effects in 
M4 is a "method of closest neighbours," where the reddening 
displacement for each program star is determined as the me- 
dian value among the closest 10 to 30 neighbours for which 
the residuals have been determined. With this approach, we 
are able to assign a reddening value on a star-by-star basis, 
allowing for any arbitrary dust geometry. The resolution is 
automatically adjusted to the number density in a given area, 
determined by the neighbour with the largest distance which 
is used. Obviously, this method smooths over, and is there- 
fore insensitive to, reddening variations on any angular scale 
smaller than 3-6 times the typical angular separation between 
stars. In general, we adopt the method described by Milone 



et al. (201 1 1 and implement only a few small changes. 

Til order to obtain the maximal discrimination between the 
effects of differential reddening and those due to photometric 
error, we use a V vs. V - K CMD, which — besides the large 
dynamic range of this colour compared to its uncertainty — 
also provides another important advantage over other filter 
combinations. In general, the direction of the reddening vec- 
tor in the CMD depends upon the value of Ry that describes 
the reddening law. This filter combination, however, is highly 
insensitive to the reddening law due to the fact that the ex- 
tinction in Ks barely changes with the dust type. Therefore 
the effect on both axes is dominated by the change in V and 
almost cancels out, with the result that we are able to perform 
the reddening correction without knowing Ry precisely at this 
point. To define the angle of the reddening vector, we use 
Ak/Av = 0.123, assuming Ry ~ 3.6, but a standard assump- 
tion of 3.1 would only change the angle of the vector by less 
then 0.5 degree. 

In addition to a reference star sample on the main sequence, 
we select a second sample along the RGB which is smaller in 
number but shows stronger sensitivity to spatial differential 
reddening because it is more nearly orthogonal to the redden- 
ing vector. 

Since our data have low photometric uncertainty only for 
specific regions in the field of M4, we exclude only stars with 
especially large uncertainties to be able to cover the majority 
of the surface area (see Table[T]). 



TABLE 1 

Selection criteria for reference stars. 



Ai'ea 


V 


ct(V - K) 


# Stars 


MS 


16.7- 18.3 


< 0.04 


4060 


RGB 


13.6- 15.7 


< 0.03 


444 
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Fig. 5. — Reddening map of M4. Blue areas indicate a lower reddening 
compared to the mean value and red areas indicate a higher reddening. Some 
differential features are presumably caused by PSF variations within observa- 
tion frames and are not due to differential reddening (see the text) The coordi- 
nate system is normalized to an arbitrary zero point at RA = 16''23™15.12'\ 
Dec = 26°25'15.4". The cluster centre is located at (x,y) a; (220, -400). 



To obtain the best angular resolution in the dense central 
regions and simultaneously be able to define diff'erential red- 
dening values for stars in sparse outer regions, we apply four 
different levels of significance: In a first step, the 30 clos- 
est neighbours to each star are selected. In case not all of 
them fall within a maximum radius of Rmax - 40", the pro- 
cess is started again, this time searching only for the closest 
25 neighbours. If necessary, the iteration is repeated for 20, 
15, and finally 10 closest neighbours until all of them are lo- 
cated within Rmax- With this approach, we assure the actual 
astrometric proximity of the closest neighbours to each of the 
selected program stars by accepting a poorer standard error 
of the mean reddening for values derived in sparse outer re- 
gions where fewer neighbours can be found in the local en- 
vironment. Any stars for which fewer than 10 neighbours are 
located within the critical radius are not assigned a differential 
reddening value. 

3.2. Results 

To visualize the spatially differential reddening across the 
cluster face of M4, we bin the surface in square cells of 
20" X 20" and calculate the median reddening residual for 
all stars that fall in this coordinate range, and for which a red- 
dening value has been assigned in the previous steps. Here, 
blue cells indicate regions with a reddening below the overall 
mean and red cells regions with values above, while the colour 
intensity increases with increasing absolute value. Cells con- 
taining only undetermined stars are shaded dark grey. Cells 
saturate when the differential reddening value becomes larger 
than |A£'(Z? - V)\ - 0.05. The resulting reddening map is 
shown in Figure |5] 

We are able to determine spatially differential effects for a 
total area of about 10' x 10' around the centre of M4. We 
estimate the size of the differential reddening by evaluating 
the median values for stars within cells of 40" x 40" and find 
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Fig. 6. — Direct comparison of our HQ-sample of M4 before (left panel) 
and after (right panel) applying the derived differential reddening correc- 
tion. Note especially, that the MSTO, the RGB, and the HB are significantly 
sharper in the coirected version of the CMD. This will help us to define the 
distance and age of M4 with a higher precision. The black bar in each panel 
indicates the direction of the reddening vector and therefore defines the di- 
rection of the shift for stars with different individual extinctions. Its length 
reflects the shift caused by a total reddening of E(B -V) = 0.37. 




Fig. 7. — Expanded view of the MSTO region before (left panel) and after 
(right panel) the differential reddening correction has been applied. We are 
able to decrease the observational scatter to about half of its original amount 
when we correct the photometry for spatial differential effects. As in the 
previous plot, the black bars indicate the direction of the reddening vector. 



a total range of k 0.2 mag between the lowest and highest 
reddening values within the area analyzed. That is more than 
half of the total mean reddening o f M4 which we determine 
to be E{B - y) = 0.37 (see section [53]l. 

Some of the spatially differential effects show strong ev- 
idence for systematic photometric errors within individual 
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frames. In particular, the vertical feature at x » in Figure l5] 
falls along the boundaries of the two most frequently observed 
regions. These systematic photometric errors are most likely 
the consequence of the PSF variations within SOFI that were 
mentioned previously. Alternatively, they could be the result 
of inadequate flat-fielding of the individual images or possibly 
crowding-dependent systematic errors in the 2MASS catalog 
itself (their angular resolution of ~ 2" is appreciably poorer 
than that of our images). 

Finally, we use the individual spatial differential reddening 
values to correct our photometry on a star-by-star basis. Fig- 
ures |6] and |7] show a comparison between our original CMDs 
and those corrected for differential reddening. With this cor- 
rection, we are able to decrease the scatter in our data by about 
50%. The reduction in the scatter along the RGB and the 
HB is especially dramatic, which allows us to determine cru- 
cial parameters like the HB luminosity or the mean value of 
E(B - V) with higher precision. 

After the correction for differential reddening, we deter- 
mine the fiducial sequence (including an empirical ZAHB se- 
quence) for our later analysis of M4. For NGC 6723, corre- 
sponding information is assessed directly from its HQ-sample 
without an additional correction for spatial reddening varia- 
tions. The fiducial sequences for both clusters can be found 
in Appendix [a] 

4. AN OBJECT-SPECIFIC REDDENING LAW 

A reddening law describes the functional dependence of 
interstellar extinction with observed wavelength. With in- 
creases in the total amount of extinction that an observed 
object suffers from, minor discrepancies in the applied cor- 
rections for reddening become more significant and can lead 
to incorrect implications for the distance or for basic atmo- 
spheric and chemical parameters derived from photometry. 
Such discrepancies can arise either from inappropriate as- 
sumptions for the actual bandpass that is used (as defined 
by the instrument and the energy distribution of the observed 
flux) or from the applied reddening law as defined by the ex- 
tinction properties of the dust in the line of sight to the object. 
In this section we will focus on the extinction properties of 
our broadband filters for a given reddening law, while we will 
address the issue of the actual shape of the reddening law it- 
self in section|5] 

The extinction along the line of sight to a source is defined 
by how much the observed flux (F^) is suppressed relative to 
what it would have been if there were no dust (F^})- The frac- 
tion of the flux which reaches Earth, passes through the filter, 
and is measured by the detector is determined by the atmo- 
spheric transmission, the mirror reflectivity of the telescope, 
the transmission of the filter and other optical elements in the 
instrument, as well as the detector quantum efficiency. The 
latter are source-independent factors and, after accounting for 
the path length through the atmosphere, can usually be con- 
sidered constant for a given observing site. Their product may 
be called system response function, which we will designate 

However, to calculate the extinction in a broadband filter 
(Aa), the energy distribution of the flux itself must also be 
taken into account. The spectral energy distribution (S a{^)) 
reaching Earth from the observed object changes with source- 
dependent factors such as temperature (Teff), surface gravity 
(log^), metallicity ([Fe/H]) and the amount of extinction (e.g. 
Ay). Knowing the system response function and the individ- 
ual stellar flux, the extinction for a given filter is given by 



AA--2.51og(^l, 



(2) 



and therefore 



-2.5 logi 



/ DAA)Ta(A)S AA)10-^'''^ ' dA 


CO I 

fDAA)T,iA)SA{A)dA ' 



(3) 



Ax describes the strength of the extinction as a continuous 
function of wavelength and — in our case — the dust-type pa- 
rameter/?!/, transforms the energy distribution of the spec- 
trum from flux units to actual photon counts in the CCD and 
therefore is equal to A in our case. Since 5^(/l) appears in both 
integrals, the physical constant involved in converting from 
energy to photon counts — he in £ = hcjA — can be extracted 
from the integral and cancel out. In the next section we will 
specify the different components of Eq. ([3]l in detail. We use 
A as a label representing some arbitrarybroad photometric 
bandpass. 

Practically, most reddening corrections are made by assum- 
ing a constant value for the extinction in broadband filters for 
every star in the sample calculated from an effective wave- 
length of the filter-detector system corresponding (in most 
cases) to a Vega-like flux distribution. However, this is only 
a first-order approximation to a more complex dependence on 
the aforementioned atmospheric and chemical parameters — 
notably the stellar surface temperature, as pointed out in for- 



mer studies (e.g. ^Bessefl et al. 1998 ; McCall 2004] [Girardi 
et al. 2008 ) — and is only valid if either the set of filters is suf- 
ficiently narrow, or the interstellar extinction is low, neither of 
which is the case for our data for M4. 

Here, we want to calculate the precise total extinction in 
each of our optical and infrared filters, appropriate to the spe- 
cific Teflf, log g and [Fe/H] for each star in a system like M4. 
The correction factors so derived can be tailored to the rel- 
atively cool, metal deficient GC stars and we want to know 
how big the difference is compared to the standard extinction 
values prov ide d in the litera ture (e.g. McCall 2004; Schlegel] 
|et al.|19 98 ; and CardelU et al.|1989| ), which are mostly derivea 
from significantly hotter 0-, B- or A-type stars. Since the 
determination of the dust type depends upon the zero points 
for the transformation factors and consequently on the stellar 
properties, these results are necessary to provide appropriate 
zero points for the reddening law of M4 in order to minimize 
systematic errors in the determination of Ry in the next sec- 
tion. To calculate the actual extinction for each filter, we use 
Eq. (|3]l and vary the shape of the targeted SED (S x) and the 
reddening law {Ax, see next section). 

One important aspect of our study will be to compare the 
colour excess in different optical and NIR filters in order to de- 
rive a consistent reddening estimate, E(B-V), for M4 from all 
filter combinations, or to explain the discrepancies. For this 
issue it is convenient to describe the reddening law in terms 
of the resultant transformation factors between the reddening 
in {B - y) and any other filter combination (Ai - A2). Those 
factors can be calculated from the extinction values for each 
filter: 



lA, 



•Aa, 



(4) 



8 



Hendricks et al. 



4.1. Ingredients for the Reddening Law Calculations 

The NIR data in J and Ks are calibrated to stars of the 
2MASS system. For this reason we use s ystem response func - 
tions for the 2MASS fihers provided in |Cohen et al.| ( |2003| l. 
These authors already combined all contributions froni the op- 
tics, detector, location-specific atmosphere and filters so that 
we can use their so-called "spectral response functions" di- 
rectly as our T,\ . 

For the optical filters and the Landolt standard stars, such a 
sophisticated synthesis of detailed characteristics is not pub- 
lished; consequently we had to calculate the response function 
using Landolt's tabulated filte r transmissions and photocath- 
ode sensitivity (Lan dolt||1992l. We adop ted the atmospheric 
transmission function from ] Allen j(| 1965, see p. 126) scaled 
to 0.7 atmospheres, approximately correct for an observatory 
2 km above sea leve l. We further included two reflections 
from aluminum from |Allen| (| 1 965| see p. 108) as is appropri- 



ate for a Cassegrain instrument. 

One of the most crucial components of Eq. (|3]l is the de- 
pendence of the interstellar extinction on the wavelength. We 
are using th e empirically derived analytical expression from 
CardelU et al. ( 1989) and specifically their equations (1), (2a), 
(2b), (3a), and (3b) to describe the extinction function. Since 
equation (1) has Ry as a free parameter, we are able to derive 
extinctions for diff'erent types of dust. 

With the use of different stellar spectra we are able to pre- 
dict the reddening law for any combination of stellar tempera- 
ture, surface gravity and metallicity, and investigate the differ- 
ential eff'ects for a stellar sample found — for example — within 
a globular or open cluster The basic problem in the selection 
of stellar spectra is the difficulty of finding a database cover- 
ing all wavelengths from near UV (U filter at ~ 3 600 A) to 
the Ks filter at ~ 21 500 A in the NIR. In total we use thi-ee 
different spectral databases in our analysis, which should im- 
prove the reliability of our predictions: 



The library of ATLAS9 model atmospheres ( Castelli & 
Kurucz 2003) 1 contains synthetic stellar fluxes for a fine 
grid in T^s and log g. Furthermore, it provides synthetic 
fluxes for stars of different metallicities from [Fe/H] - 
+Q2 down to [Fe/H] = -4.0. 



The atlas from |Pickles| ( | 1998) 1 contains observed stel- 
lar fluxes, compiled by combining the data from sev- 
eral existing catalogues overlapping in wavelength cov- 
erage. The atlas provides 65 flux-calibrated spectra 
spanning a wavelength range of 1 150 - 25 000 A and 
therefore encompasses all relevant optical and NIR fil- 
ters. Since we expect the reddening correction to de- 
pend on both dimensions (i.e., Teff and log g) of a CMD, 
we select a sample of spectra which follow a typical 
isochrone representation of a GC fiducial sequence as 
closely as possible. The vast majority of the spectra 
have solar metallicity ([Fe/H] = 0.0) so that we are not 
able to make predictions for metal deficient stars using 
this atlas. 



The MILE S Library of Empirical Spectra ( Sanchez 
Blazquez et al. 2006) provides observed spectra for 
stars with a wide range of metallicities and atmospheric 
parameters. We were unfortunately not able to select 
a sufficiently large sample of equal-abundance stars to 
occupy the whole evolutionary sequence of a GC. Fur- 
thermore, the wavelength coverage ranges from 3 525 A 



to 7 500 A and therefore fully covers only the B and 
the V filters. For this reason the MILES library is used 
only to empirically test the results for low metallicities 
as derived from synthetic spectra (see section [4.2. 3[ ). 
For this purpose we select a sample of spectra located 
close to the MSTO location in a CMD (T^ff k 6000K; 
\ogg X 4.5) with different metallicities typical of GCs. 

4.2. Results 

With the extinction law defined by Eq. (|3]l and the different 
components described in section |4 1| we mvestigate the im- 
pact of different GC parameters on the reddening law which 
has to be applied to correct the data for interstellar reddening. 
In most cases we discuss the results using as an example the 
transformation factor Fy^K between the observed reddening 
in B - y and the derived colour excess inV - K. In this way, 
the impact on a CMD that combines optical and near infrared 
filters can be directly assessed, whereas a discussion on the 
basis of Aa would keep the reader in the dark about the ef- 
fects found when several extinctions are combined according 
to Eq. Q. Throughout this section we assume a standard dust- 
type parameter of Ry - 3.1. A higher value of Ry, however, 
would only cause a constant shift in F^^^a^ but does not affect 
the relative changes caused by variations in stellar parameters. 

Throughout this paper we will use the expression 
reddening-law zero point for the value of Fa, -a, or A a that 
is obtained for a star in the TO region of the CMD. It is there- 
fore equated to the result of Eq. ([3]l for a star with atmospheric 
parameters Teif - 6000K and \ogg - 4.5. This definition has 
been adopted since, in most cases, this part of a GC evolu- 
tionary sequence is used to fit isochrone models to the data. 
In the case where the zero points are directly referred to M4, 
they further imply a metallicity of [Fe/H] = -1.0. In contrast 
to that, we discuss the differential effect in /^Ai-A2 towards a 
given zero point due to variations in stellar parameters or — in 
the following section — Ry. 

4.2.1. Effects of Temperature 

Figure |8] shows the effect of varying the temperature on the 
transformation factor Fy^j^. The results obtained with ob- 
served spectra are shown in direct comparison to those from 
synthetic spectra. For this comparison we choose a stellar 
sample from the Pickles atlas which approximately follows 
the MS of a GC evolutionary sequence, while for the ATLAS9 
synthetic spectra we choose a constant log g - 4.5 and a wide 
range of temperature. 

Since we observe no significant differences in shape for the 
filter combinations relevant for this work (including B, V, I, 
J, H, or Ks), the most significant features of Figure [8] can be 
generalized to: 

1) The results obtained with the observed spectra show 
a relatively smooth increase for Fai-a^ with decreas- 
ing temperature. Between the TO temperature and the 
coolest stars seen on the lower MS at log T^s ~ 3.58, 
the change in Fa,-Ao is about 3%. 

2) The synthetic spectra match the results from the 
observed sample very well for temperatures above 
5 000 K, whereas there is a strong deviation for lower 
temperatures. For T^ff < 5 000K the synthetic spec- 
tra display a strong decrease in the transformation fac- 
tor after reaching a maximum value, whereas the ob- 
served spectra show a monotonic increase of Fai-A2- 
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Fig. 8. — Dependence of tlie transformation factor Fy-K on stellar temper- 
ature. In the upper panel, diamonds indicate spectra from the observed atlas 
from Pickles 1 1998i and black li nes show the res ults obtained with synthetic 
ATLAS9 spectra from Castelli & Kurucz| j 2003) . The grey area indicate the 
range of typical GC temperatures with the lower boundary at temperatures 
for cool MS stars and the upper boundary at roughly TO temperatures for 
old stellar populations. The lower panel shows the differences between the 
models and real spectra at the solar metallicity. Good agreement is found for 
temperatures above ~ 5 000 K. 



This discrepancy seems to be stronger for more metal 
rich stars since the synthetic spectra with a metalhcity 
of [Fe/H] = -2.0 yield almost the same shape as the 
observed spectra. We further investigate the discrepan- 
cies between observed and synthetic spectra in the next 
section and discuss possible reasons for them. 

3) Within temperatures ranging from 4 500 K to 6 500 K, a 
decrease in metallicity causes a significant decrease in 
the zero point of Fa, -a, whereas the shape or the differ- 
ential effect is almost the same. For Teff > 9 000K al- 
most identical values of /^Ai-At are found for all metal- 
licities: metallic features have no perceptible eff'ect on 
stellar spectral energy distributions at these tempera- 
tures. 

4) Principally, our zero points for the transformation fac- 
tors which include J, H, or Ks filters seem to be higher 
then the ones commonly used in the literature. For ex- 
ample, Fv-K = 2.78 - 2.85, depending on the metallic- 
ity, whereas t he canonical v alue using a Vega-like star 
is about 2.71 ( |McCall|2004l l. 



4.2.2. Discrepancies between observed and synthetic spectra 
According to Eq. ((4|), Fa,_a, includes the extinction for at 



least three different filters. To find the source of the discrep- 
ancy between the models and the real spectra in Figure [8] we 
looked into each component separately and compared the ex- 
tinction in all filters for the observed and the synthetic spectra 
to determine whether there are general differences or whether 
they are related to specific filters. The results are shown in 
Figure |9] We find that there is very good agreement for each 
filter except U and, more dramatically, B which explains very 




Fig. 9. — Relative extinct ions for all filters as derived from the observed 
spectra fromlPicklesV1998 l (diamonds) and synthetic ATLAS9 spectra (solid 
line) from Castelli & Kurucz (2003 1: Although the actual numbers are dif- 
ferent, the y-axes show the same scale for each subplot (Aa/Av = 0.0075 
between tickmarks) so that relative changes with temperature and dift'erences 
between synthetic and observed spectra can be compared directly. There is 
good agreement for all photometric filters, including the absolute value for 
the V-band. Only the extinction in the B-filter shows significant differences 
between the different types of spectra for temperatures below 5 000 K. The 
range of typical GC temperatures lie between the dashed vertical lines. 



well the deviation between the Pickles and ATLAS9 spectra 
in Figure |8] because a higher value of results in a lower 
value of Fy-K- Since Ab contributes to every transformation 
factor, it can explain the observed deviations for all filter com- 
binations. 

To explain the difference in the spectral energy distribu- 
tions, we overplotted a real spectrum with the synthetic spec- 
trum in the wavelength range of the B filter for T^s - 3 500 K, 
where we expect a strong deviation, and for T^s - 5 000 K 
where we expect little or no deviation (see Figure [T0| . In 
the cooler star, the synthetic atmosphere models (blackcurve) 
consistently overestimate the true stellar flux (grey curve) by 
many percent at the shorter wavelengths within the bandpass, 
while they consistently Mniierestimate the flux in the longward 
portion of the bandpass. This causes a shift in the photon- 
weighted mean wavelength of the light that passes through 
the filter (vertical lines), which results in a mis-estimate of 
the net extinction integrated over the B bandpass. 

For low temperatures, molecular lines and bands become 
more prominent in the optical part of the spectrum, includ- 
ing the wavelength range of the B filter Specifically, the 
wavelength range encompassing the observed discrepancies 
includes prominent CN and CH molecular bands which are 
widely used as spectral indices to investigate variations in 



the CNO-cycle elements in GCs (Norris & Freeman 1979[ ). 
The indices S3839, S4142 ( for CN ) and CH 4300 (for CH) 
as defined in Harbeck et al. (2003) are located exactly in the 
region where we ob serve the strongest deviations (see e.g., 
Pancino et al.|2010 Fig. 3). Apparently, the synthetic atmo- 
sphere mo dels do not per fectly reproduce these complex fea- 
tures (e.g., Bessell et al.| [r998 ). This fact, together with the 
increasing deviation towards lower temperatures and higher 
[Fe/H] supports the assumption that the results obtained with 
the synthetic ATLAS9 spectra are not entirely trustworthy for 
the cool temperature regime due to an inadequate treatment of 
molecular bands. 
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Fig. 10. — Detailed plot of the different spectral energy distributions at 
Teff = 3 500 K (top) and T^s = 5 000 K (bottom) to visualize the cause of 
the extinction differences in the B-band. Vertical lines indicate the resultant 
effective wavelength for B when the different spectra are used (thick black: 
synthetic ATLAS9; grey: observed spectrum). The significantly shorter ef- 
fective wavelength for ATLAS9 spectra in the top panel results in the higher 
relative extinction in B compared to observed spectra. For T^g = 5 000 K. 
The differences in the spectral flux are much smaller and the effective wave- 
length is basically the same for synthetic and real spectra. Both spectra are 
normalized to the same wavelength so that intensity differences can be com- 
pared directly. 



4.2.3. Effects of Metallicity 

Since the different spectral samples agree very well for 
temperatures above 5 000 K, we infer that the synthetic mod- 
els predict the change of /^a.-a, zero points with metallic- 
ity correctly at MSTO temperatures . To confirm this infer- 
ence, we use the spectral atlas from'Sanchez-Blazquez et al. 
(|2006) to compare the results of the ATLAS9 models with 
the observed spectra for a wide range of metallicity between 
[Fe/H] = -1-0.2 and [Fe/H] = -2.5 at a constant temperature 
of Teff ^ 6 000K. 

Unfortunately, we are not able to generate extinctions for 
each filter with the observed metal poor spectra, due to their 
limited wavelength coverage. For this reason we derive trans- 
formation factors where the observed metal poor spectra are 
used only for the B and V filters, whereas we derive the ex- 
tinction for the remaining bands with synthetic metal poor 
spectra. Then, we compare this hybrid model to a purely syn- 
thetic one where the extinctions for all filters are derived from 
synthetic spectra. Note that for Fy-x, for example, only the 
extinction in the Ks filter is calculated with synthetic fluxes. 
Since we expect changes in the metallicity will mainly im- 
pact the optical part of the spectrum (and therefore U, B, and 
V), we are confident that the results obtained with the hybrid 
model reflect the true situation very well. 

Table l2] summarizes the parameters of Sanchez-Blazquez 
and the ATLAS 9 sp ectra which are used for the direct com- 
parison. Figure 1 1 shows the derived values for Fy-x using 



the synthetic anothe hybrid model, together with the dif- 
ferences between them. The consistency between synthetic 
metal poor spectra and the observed ones is extremely good. 
Only the results from the most metal rich synthetic spectra for 



[Fe/H] > 0.0 deviate from the values obtained with observed 
fluxes. For the rest of the direct comparisons the difference in 
Fv-K is considerably smaller than the photometric uncertainty 
without any systematic deviation. 

Reviewing Figure [8] we want to point out that the solar- 
metallicity synthetic spectra reproduce very well the shape of 
all fA,-A2 derived from observed spectra above 5 000K. By 
comparing synthetic spectra with different metallicities, we 
find the shape to be very similar for all metallicities in the 
temperature range 5 000 K < T^ff < 9 000 K. Therefore, a 
change in metallicity in this temperature range mainly pro- 
duces a shift in the zero-point value of the transformation fac- 
tor. As a result, we can calculate these zero points from metal 
poor synthetic spectra in a temperature range where the re- 
sults are confirmed by metal poor observed spectra and then 
apply them to different temperatures. However, it is important 
to point out that the effect of metallicity on the reddening law 
is not constant over larger temperature scales and the use of 
a constant zero-point offset is an approximation that seems to 
be justified only for the limited range in GC temperatures. 



TABLE 2 

Parameters for observed and synthetic spectra in direct comparison. 
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Fig. 11. — Top: Zero points for Fy-K as a function of [Fe/H] for the hy- 
brid combination (diamonds) and a pure synthetic sample (circles) are com- 
pared directly. The dashed line indicates tlie functional dependence of Fy-K 
with [Fe/H] as derived from pure synthetic spectra with Teff = 6 000 K and 
logg = 4.5 which represents the atmospheric parameters at the MSTO. Note 
that these parameters are not used for all ATLAS9 spectra in the direct com- 
parison sample (see Table [2| and therefore some circles do not fall exactly 
on the dashed line. Bottom: Deviation in Fy-K between pure synthetic and 
hybrid metal poor spectra. AFv-a: = Fy-KUyn) - Fv-K^hy^)- 
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Fig. 1 2. — Effects of surface gravity on the transformation factor Fy-K for 
a wide range of temperatures. For tliis calculation we used synthetic spectra 
with [Fe/H] = -1.0 and [a/Fe]= +0.4 to simulate the chemical composition 
of the stars in M4. As in previous plots, the grey area indicates the range of 
typical GC temperatures. In general, a lower surface gravity results in higher 
values for Fy-K, where this effect seem to increase with lower temperatures 
and is not observable for temperatures near the MSTO. The effect seems to 
be inverted for temperatures hotter than T^g » 6 300 K. 



4.2.4. Effects of Gravity 

The absolute magnitude My of a star in a simple stellar 
population such as a GC system is strongly correlated with 
its surface gravity and therefore has a possible impact on the 
stellar spectrum and the reddening law which has to be ap- 
plied. Figure 12 shows the effect of different \ogg values 




on the transformation factors derived with synthetic spectra 
from ATLAS9. Ignoring the temperature range cooler than 
the peak in Fy-K, we observe an increase for /^a,-a, with de- 
creasing surface gravity. The effect becomes more significant 
towards cooler temperatures, where it is of the same order as 
the temperature effect itself. For the lowest GC temperatures, 
the change in Fa, -a, between a dwarf and a giant is about 
3% for all transformation factors. This is especially important 
for stars on the RGB since the decrease in temperature and 
surface gravity combine and enhance the differential effect of 
reddening along the RGB sequence. However, for tempera- 
tures between 5 600 K and 6 300 K the differential effect with 
gravity practically disappears and eventually reverses for even 
higher temperatures. 

4.2.5. Effects of Extinction 

The total strength of interstellar extinction itself is a factor 
which impacts on the effective wavelength of the filters. This 
is because interstellar extinction is a function of wavelength 
and therefore changes the distribution of the incoming flux by 
damping the blue part of each bandpass more strongly than the 
red part. Figure[T3]shows the temperature-dependent value of 
Fv-K for different absolute extinctions Ay- Compared to at- 
mospheric parameters, the influence of A v on the transforma- 
tion factor is very small and although we use E{B -V) - 0.36 
in our calculations to tailor the results for M4, they are basi- 
cally valid for at least all E(B -V)< 0.5. 
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Fig. 13. — Effects of total extinction for four different total extinctions cor- 
responding to E{B-V) = 0.1, 0.3, 0.5 and 1.0 assuming Ry = 3.1. The effect 
is small compared to atmospheric influences that have been discussed in the 
previous sections and has about the same impact for all temperatures. The 
grey area indicates the range of typical GC temperatures. 



4.2.6. Star-by-Star Variations within a Stellar Population 

Despite the definition of appropriate zero-point values for 
relative extinctions Aa/Ai/ and transformation factors Fa, -a, 
appropriate to the atmospheric and chemical conditions of GC 
stars, the question comes up whether changes of those pa- 
rameters within a typical GC population cause any significant 
change in the appearance of the cluster fiducial sequence com- 
pared to a constant correction with respect to the photometric 
uncertainty. Such changes could affect the determination of 
[Fe/H], distance, age and other parameters when determined 
from photometric studies. 

For that reason we compared the extinction for a faint MS 
star (Feff ~ 3 680 K, My ~ 9.05) to a TO star (T^ff ^ 5 810K, 
My ^ 4.25) and a stai- on the upper RGB (Fefj ~ 3 820K, 
My = -0.66) to simulate the highest possible differential ef- 
fects in both temperature and surface gravity for a reddening 
of E(B-V) = 0.36. Using SEDs with appropriate parameters, 
we found differential effects of AE{V - K) ^ 0.025 mag be- 
tween the MS stai- and the TO stai- and AE(V -K) ^ 0.05 mag 
between the RGB star and the TO star. This effect is small 
and, for both examples, of the same order as or even below 
the photometric uncertainty. Therefore the shape of the evo- 
lutionary sequence of clusters with significant reddening such 
as M4 are not affected significantly by differential effects due 
to variations in temperature and surface gravity. That means 
that parameters like metallicity and age — when derived by 
isochrone fitting — are not significantly biased by reddening 
zero point variations within the population. Specifically, the 
luminosity of the HB, which is often used to derive the dis- 
tance modulus of GCs, is not biased when the zero-points 
are normalized to TO temperatures; at those temperatures, the 
change in surface gravity does not affect the reddening zero 
points. In general, however, the error caused by the use of 
a constant reddening correction is of systematic nature and a 
star-by-star correction should be considered when the highest 
possible photometric precision is desired or if the reddening 
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in the line of sight is larger than for M4. In Appendix [A] we 
therefore provide two fiducial sequences of M4: one far the 
(apparent) observed sequence and one where the aforemen- 
tioned differential effects are corrected on a star-by-star basis 
with respect to MSTO conditions. 

4.2.7. Extinction Zero Points 

From our calculations, we provide zero points for relative 
extinction values (A a /Ay) for a wide range of different tem- 
peratures, metallicities, and values of Ry- From them, ap- 
propriate colour transformation factors can be derived. The 
metallicity and Ry tables have been created with synthetic AT- 
LAS9 spectra and observed fluxes are used to calculate val- 
ues for the temperature table due to the uncertain treatment 
of molecular absorption in the models. For all tables, we as- 
sume Ay - 1.12 (suitable for M4), but the values are valid for 
practically all Av < 2. These transformation factors should 
be applied to correct the colour excesses in filter combina- 
tions of optical Johnson-Cousins and/or NIR 2MASS filters. 
Since the tabulated differential effects due to temperature and 
[Fe/H] are basically independent of Ry, the tabulated offset 
due to the dust type can be used for any combination of the 
atmospheric and chemical parameters. The same can be as- 
sumed for the Tsff and [Fe/H] effects in the range between 
5 000 K and 9 000 K. 

The extinction tables are listed in Appendix [B| where they 
are compared to pr ominent l iteratu re values. Note that for 
J, H, and K„ only |McCall| p004| uses the actual 2MASS 
bandpasses, whereas all other references use slightly different 
filters. 

4.3. Comparison to existing Work 

Some studies have considered the effects of temperature, 
surface gravity or other parameters on the reddening law. 
However, most of them use onl y synthetic spectra. The recent 
work from llGirardi et al.| ( |2008] l (hereafter G08) uses ATLAS9 
spectra (the same generation than we adopt) to produce a dif- 
ferential reddening law with a grid of log g and Teft specif- 
ically applicable for the WFPC2 and ACS photometric sys- 
tems from the Hubble Space Telescope. They found a differ- 
ential effect of tricolour) ~ 0.05 mag for a total hypothetical 
extinction of Ay x 3.0 mag, which corresponds to a redden- 
ing of E(B - V) ~ 1.0 using a standard Ry value of 3.1. We 
find a similar sized effect for less than half the reddening us- 
ing aV-K magnitude plane, which can be explained with the 
larger wavelength range that is spanned by this filter combi- 
nation. All values from G08 are derived with solar metallicity 
spectra only; they claim that the variation with metallicity is 
negligible. That is not what we find in our study: even if 
the metallicity effect for the absolute extinction is very small 
in a single filter, it can increase significantly when combining 
four filters to derive the transformation factors between differ- 
ent colours. Using the same thresholds for [Fe/H] as G08, we 
find a difference of SFy^x » 0.08 mag between [Fe/H] = 0.0 
and [Fe/H] - -2.5, a change of 2.8% which is about ten 
times larger than G08 states. Overall, our work reveals a 
roughly equally sized effect for changes in [Fe/H], T^ff, and 
\ogg {~ 3% each) encompassed by GC stars, whereas G08 
claim that the effects of temperature variations dominate over 
those of surface gravity while those arising from metallicity 
variations are negligible for the HST filters. 

Another interesting issue to discuss is the quality of AT- 
LAS9 spectra, especially for low-temperature stars where 



molecular lines become important. Our work finds very good 
agreement between observed and synthetic s pectra for aU 
bands except B for Teft < 5 000 K. Interestingly, Bessell et al. 
( 1998) use an older version of ATLAS9 models (from 1993) 
and found a good fit only for temperatures above 4 500 K due 
to their incomplete molecular opacity at lower temperatures. 
Although we are using a more recent version (2003) with an 
improved treatment of molecules, we seem to find a similar 
effect setting in at almost the same temperature. Carefully ex- 
amining Fig. 2 in G08 indicates that their study suffers from 
the same problem: the relative extinction A/^/Ay (As JAy in 
G08) for the F435W filter — which is closest to Johnson B — 
shows a dip at cooler temperatures similar to our results with 
ATLAS9, whereas the use of observed spectra reveals a mono- 
tonic decrease in relative extinction towards lower tempera- 
tures. 



The work from McCall| ( |2004| l is primarily focused on 



the extinction correction of mtegrated SEDs for extragalac- 
tic sources, suitable — for example — for the determination of 
accurate distances from Cepheids. However, the author dis- 
cusses the effec t of temperatur e for resolved stellar spectra as 
well, using the Pickles ( |1998| l dataset and found differences 
up to 23% between unevolved 05 and M6 stars at the so- 
lar metallicity and a nearly equal effect between very cool 
evolved and unevolved stars, which is consistent with what 
we find when we consider the full range of temperature. Note 
that the largest effect is caused only by the coolest stars with 
Teft < 3 500 K and therefore differences are much smaller if 
limited to the range of GC temperatures. Unfortunately, Mc- 
Call only provides zero-point values normalized to a Vega- 
like flux together with Ry - 3.07 for a "normal" reddening 
law, which makes it difficult to apply his results to globular 
cluster conditions and compare them to our results in more 
detail. 

5. THE DUST TYPE OF M4 

The type of dust that causes the interstellar reddening can 
show significant variations for different lines of sight and is 
responsible for significant changes in the properties of ab- 
sorption and scattering. Therefore the actual reddening law is 
mainly determined by the specific properties of the dust. This 
makes it, in principal, necessary to determine an individual 
reddening law for each object that suffers significant interstel- 
lar extinction, before correcting for those effects appropriately 
without inducing systematic errors. To good approximation, 
the dependence on the dust type can be parameterized by the 
ratio between total-to-selective extinction, Ry = Ay/E(B-V), 
which seems to be strongly correlated with the gr ain proper- 
ties of the dust ( jCardem et al.|1988|pathis|1990| l and varies 
between Ry 2.6 and a value up to ~ 6.0 for dense, cold 
molecular clouds or star forming regions. Since it defines the 
shape of the reddening law and the ratio between extinction 
and reddening, a precise knowledge of Ry (or some equiva- 
lent parameter) is crucial for the absolute distance determina- 
tion of stellar systems. Usually a value of ~ 3.1 is used to 
describe the foreground reddening for GCs in the Milky Way 
as it is associated with the average extinction properties of the 
ISM. However, there is strong evidence in the literature that 
M4 follows a different, anomalous reddening law due to the 
fact that the cluster lies behind the Sco-Oph dark cloud com- 
plex which adds an additional uncertainty to the distance and 
possibly to the age of the cluster as well as to the assumed 
mean reddening expressed in E{B - V). 

In this section, we want to determine the properties of the 



A New Reddening Law for M4 



13 



dust in the line of sight to M4, characterized b y the appropri- 
ate val ue of Rv in the reddening law given in Cardelli et al. 
1989[ ). From this, the actual ratio between Ay and E{B - V) 



for our photometry can be determined in order to derive the 
distance to the system. 

5.1. Method 

In general, we want to find out which value of the dust type 
parameter Ry produces a consistent match for isochrones with 
the observed CMDs of M4 in all filter combinations once the 
empirical offset is given for one colour Here, we make use of 
the possibility of treating the dust type as a free input param- 
eter to calculate transformation factors Fa, -a-, between dif- 
ferent colours as a function of Ry- The dust type parameter 
enters equation ([3]) twice: First, it defines the slope of the 
reddening law according to the Cardelli et al. equations. Sec- 
ond, it defines the actual extinction Ay for a given reddening 
E(B - V). 

As far as we are aware, the only similar attempt to deter- 
mine the ratio of total to selective extinction in the line of sight 
to star clusters can be found in An et al., ( ,2007| ). Here, Ry is 
determined by a minimization of the shape of the redden- 
ing law defined by the observed colour excesses in optical and 
NIR broadband filters. In contrast, we choose one reference 
colour (B - V) and find the best value of Ry by matching all 
other filter combinations individually. Compared to An et al. 
(|2007 ), our approach is more sensitive to variations in the as- 
sumed value of [Fe/H], mainly because metallicity changes 
the theoretical colours for the isochrones and therefore the 
adopted value of E{B - V). On the other hand, we are able to 
minimize filter-specific systematic errors for our bandpasses, 
while a minimization of the overall shape does not take 
those distorting effects into account and therefore can lead to 
systematically wrong answers if in fact the theoretical colour 
for only one bandpass is inaccurate. Since we have two pow- 
erful tools to constrain the metallicity of M4 (spectroscopic 
data and CMD isochrone fitting) but almost no access to the 
accuracy of synthetic filter colours, we decided not to use a 
minimization algorithm to determine Ry. Furthermore, An 
et al. indeed account for temperature effects in the broadband 
filter extinct ion properties by us ing the colour-correction term 
provided by Bessell et al. ( 1998 ), but they do not consider any 



adjustment for gravity or metallicity as we do in our study. 

5.1.1. Empirical Color Offsets 

To determine the empirical colour offset between the ob- 
served CMD and the theoretical isochrone for each combi- 
nation of filters, we use an age-independent fiducial point on 
the MS, defined to be 5.5 mag below the red ZAHB V-band 
magnitude. For M4, this occurs at V = 18.96. Our fiducial 
point is insensitive to age, since it lies well below the MSTO 
and it is further defined with respect to the ZAHB luminos- 
ity, a standard c andle which is know n to be nearly indepen- 
dent of age (e.g.. Stetson et al.|1996) l. This is important since 
ages are not precisely known for GCs and could consequently 
cause a systematic error in the determination of Ry. We fur- 
ther choose the fiducial point intentionally below the MSTO 
to avoid known uncertainties within theoretical isochrones for 
evolutionary stages later than the TO due to the onset of — yet 
poorly understood — deep mixing processes (see, e.g., Van- 
|denBerg 2005 ). 

We find the colour of the fiducial point in the observed 
CMD by calculating fiducial sequences along the MS of M4 



in each desired filter combination using data that have been 
corrected for spatial differential reddening. Fiducial points 
are obtained by binning the data in magnitude intervals of 0.3 
along the MS and calculating the median value for each bin. 
The fiducial sequence is then defined by cubic spline inter- 
polation between the fiducial points. This approach has the 
advantage of being consistent for all the colours considered: 
since both the reference luminosity and the colour of the fidu- 
cial point are derived the same way for all filter combinations, 
any systematic uncertainty cancels out to first order. In other 
words, even if the derived fiducial sequence might not rep- 
resent the actual median reddening value for a given colour 
(e.g., due to a bias caused by binary stars), it still does rep- 
resent the same reddening value for all colours. Furthermore, 
using the median value to define the fiducial points, we mini- 
mize the influence of outliers such as field stars, blended stars 
and specifically binary stars that do not fall on the actual MS. 

Finally, each colour offset is given by the colour difference 
between the empirical fiducial point and its theoretical equiv- 
alent on the appropriate isochrone. 

From the individual observed colour excesses we derive the 
necessary transformation factors Fa, -a, for each filter com- 
bination and determine the corresponding Ry values from our 
reddening law calculations. The zero points of the transfor- 
mation factors are calculated for [Fe/H] =-1.0 (see section 
5.2 for a discussion about the metallicity of M4), E(B -V)- 
0.37 and for the atmospheric parameters Tefj = 5 250 K and 
log g - 4.5 based on the location in the CMD of the fiducial 
point that is used to determine the offset. 

5.1.2. Error Estimation 

In addition to the best fitting values of Ry, we estimate for 
each colour the observational standard error of the derived off- 
set to the isochrone {crihcolour)) and the uncertainty in Ry 
(cr(Ry)) which defines the overall precision of our method. 
Here, we have to take into account the photometric uncer- 
tainty of both the observed ZAHB luminosity (o-CVhb)) and 
the colour of the fiducial point at that luminosity (cricolour)): 



d 

(p- {decolour) - (P-{Vhb){ + cP' {colour). (5) 



To determine the first component on the right hand side of 
Eq. ([5]l, we assume criyHB) - 0.03 mag and calculate the 
change in colour of the fiducial point when its luminosity is 
shifted by that amount. For the latter term we calculate the 
median cr(colour') at the fiducial point luminosity. The quan- 
tity a-{Kcolour) thus depends mainly on the photometric qual- 
ity in each filter and to some extent on the slope of the MS at 
the fiducial luminosity. The uncertainty of Ry for each filter 
combination is then determined by the propagation of those 
errors from {(T{l^coloury) to a-(Ry). Here, the governing crite- 
rion is the sensitivity of each filter combination to Ry, which 
is discussed in the next section. Finally, each filter combina- 
tion yields an individual best fitting value for Ry together with 
an interval where the colour excess is reproduced to within its 
observational uncertainty limits. In addition to the observed 
random errors discussed here, the main sources of systematic 
errors are discussed in section l5?4l 



5.1.3. Independent Filter Combinations 

Not all filter combinations have the same uncertainty and 
not all of them yield independent results. We therefore have 
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Fig. 14. — Reddening law for different values of Ry. The effective wave- 
length of our filter set is indicated with circles and triangles for Ry = 3.1 and 
R\i = 4.1 respectively. The change of the reddening law with the dust type 
parameter is the crucial key to determine the dust type relevant to M4. 



to decide which colours we want to use to minimize the final 
value of o-(Rv)- 

The first decision which has to be made is the reference 
colour to which other combinations shall be compared, where 
the crucial criterion is to maximize the sensitivity of the ratio 
^(Ai - A2) I E{Aref\ — J^refi) to Ry- This scnsitivity is de- 
fined by the difference in the effect for the correction factor 
of a given filter combination compared to the correction fac- 
tor of the reference colour when Ry is changed. Figure [14] 
illustrates the effect of Ry on our filters. It shows the relative 
extinction as a function of wavelength, assuming two differ- 
ent values of Ry- Regarding this figure, the most promising 
reference colour is E(B - V), because its change with Ry dif- 
fers the most from that of all other combinations. From this 
perspective, E{U - V) would be an even better choice. How- 
ever, with an effective wavelength of ~ 3 560 A, the U filter is 
located on the extreme short-wavelength end both of common 
detector sensitivities and of the atmospheric transmission. It 
is therefore very sensitive to the exact atmospheric conditions 
and the instrument response. Moreover, the bandpass falls on 
the Balmer convergence and jump in stellar spectra, which 
makes the filter sensitive to even small variations in Teif and 
log g. Consequently, we expect significant uncertainties in the 
synthetic U magnitudes which would cause a systematic error 
in the determination of Ry. 

The uncertainty of every individual comparison colour de- 
pends mainly on two factors: the photometric quality of 
the data in the given filters and the sensitivity of the ratio 
E(Ai - A2)/E(B - V) with Ry. The latter factor is optimized 
by the filter combination spanning the highest wavelength 
range and the biggest difference from E(B - V) regarding Ry- 
sensitivity. Given E{B - V) as the reference value, we can 
assign a maximum of four independent colours contributing 
to the determination of the dust type under the assumption 
that results which include the U filter are not trustworthy and 
that R can be reliably predicted from V and /, so it does not 
really provide independent information. 

To minimize systematic errors within the individual band- 
passes, we choose the colours in a way such that each fil- 
ter appears on alternating sides. Therefore we use E(V - I), 
E(I - J), E{J - K), and E(B - K) as independent colour mea- 
surements to be compared to E(B - V). This choice might 
seem surprising at first, since those filter combinations all 



have relatively large individual cr{Ry) compared to E(V - J) 
or E{V - K) for example (see Table However, by using 
alternating combinations we are reducing any systematic er- 
rors from faulty colour-temperature relations (CT-relations) 
used by the isochrones or bad zero-point calibrations of the 
observed data to the standard system; the resultant high accu- 
racy outweighs the high precision of some individual colours. 
The method therefore necessarily relies on a combination of 
optical and NIR photometry, since otherwise only E{V - I) 
could confidently be used to determine the dust type, which 
would dramatically increase both the systematic and statisti- 
cal errors in the result. 

5.2. Isochrones 

The method we describe here to determine the dust type 
toward a stellar population relies fundamentally on the as- 
sumption that theoretical isochrones fit observational data ei- 
ther if the correct reddening law is used, or — of course — when 
there is no reddening. The main source of uncertainty in the 
fit of isochrone models to observed fiducial sequences comes 
from CT-relations, which predict a synthetic colour from the- 
oretical surface temperatures and gravities (from model at- 
mospheres) as well as from the predicted Teji-scales, while 
the theoretical lumi nosities are considered to be more robust 
( |VandenBerg"|2005) . 

We use me newest generation of Victoria-Regina 
isochrones to determine the colour offsets in different 
filter combinations. These models have been significantly 
updated since the latest description in VandenBerg et al. 
(12006). A summa ry of recent changes is provided by |Bro- 
[gaard et aL (201 1 1, but the most importa nt is the treatment of 
the gr avitational settling of helium (see Proffitt & Michaud 
1991| l. Besides theoretical sequences reaching from the 
MS up to the RGB tip, we use ZAHB models to determine 
{m - M)y, which is necessary to match our fiducial point 
luminosity between observed CMD and isochrones. Once the 
reddening and the dust type are known, we can further use 
the HB as a standard candle to determine the distance to the 
cluster 

CT-relations and the overall fit of Victoria-Regina 
isochrones have been tested intensively for both optical and 
near-infrared filters in tw o recent studies: for optical bands 
VandenBerg et al. ( 2010| l found very good agreement between 
Victoria-Regina isochrones and data for local Population II 
subdwarfs with Hipp arcos -based distances. In a further test, 
the isochrones reproduced the observed fiducial sequences 
for several globular clusters with different metallicities very 
well. For their study, they tested th ree different CT-relations : 
the empirical transformations by jCasagrand e et al.| (20101, 
those derived from M ARCS model a tmospher es an d the semi- 
empirical relations by Vandenberg & Clem ( 2003V for which 



they found only minor ditt'erences. According to this paper, 
the models fail to match the observed sequences only for faint 
MS stars with My > 6.5. 



Brasseur et al. (2010 1 tested Victoria-Regina isochrones us- 
ing MARCS model atmosphere-based CT-relations for the 
2MASS J, H and bands against the same sample of field 
subdwarfs and found good agreement between the models and 
the data. Only a small offset was detected in J - K and 
V - K, where the models seem to predict those colours too 
red by ~ 0.03 mag each. Comparing isochrone sequences 
to the observed CMDs of globular clusters generally yields 
good agreement for all except the lowest metallicity systems 
([Fe/H] < -2.0). Those discrepancies, however, are limited 
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Fig. 15.— Direct overlay of M4 (black) and NGC 6723 (grey). The data 
for NGC 6723 have been shifted by A(V - K) = +1.07 and AV = -2.06 to 
match the sequence of M4 (NGC 6121). Although the MS for NGC 6723 
is too poorly defined to compare the sequences directly at large magnitudes, 
both clusters show a remarkable consistency in their RGBs, SGBs, and HE 
sequences, implying that both systems are almost identical in metallicity and 
age. 



to the lower RGB segments, where the isochrones seem to be 
significantly too red, whereas the location of the MS matches 
for all metallicities. 

With the point on the MS 5.5 mag fainter than the ZAHB 
and a metallicity in M4 of [Fe/H] = -1.0, we satisfy both 
critical conditions for magnitude and metallicity very well, so 
that we are confident that the models should match the data if 
the correct reddening law is adopted. 

5.2.1. The Guinea Pig: NGC 6723 

We use optical and NIR photometry for the nearly unred- 
de ned globular cluster NGC 6723 {E{B- V) ~ 0.05 according 
Harris|[T996| l to test the morphological fit and the consis- 



to 

tency of our isochrones in different observed colour combina 
tions under optimal conditions. NGC 6723 is almost identical 
to M4 in chemistry and age and is therefore perfectly suited 
as effectively unreddened reference cluster. In Figure 15 we 
correct NGC 6723 to the distance modulus and reddenmg of 
M4 to obtain a direct overlay of the two systems. Both se- 
quences show exactly the same slope along the RGB and con- 
tinuing to below the TO, which implies almost identical ages 
and chemical abundances. On the lower MS, the photometric 
quality of NGC 6723 is too poor to allow a precise compar- 
ison. Remarkably, even the HB sequences are superimposed 
precisely with a similar bimodal distribution of red and blue 
HB stars — which is unusual, if one considers their relatively 



Fig. 16.— Fits of 11 and 12 Gyr isochrones for Y = 0.25, [Fe/H] = -1.0, 
and [ci/Fe] = +0.4 to various filter combinations of NGC 6723. Using a 
reddening of E{B - V) = 0.05 and a standard reddening law with Ry = 3.1, 
the isochrones show a good match to the observations for each combination. 



high metallicities. 

For NGC 6723, a metalhcity of [Fe/H] = -1.0, an age of 
11-12 Gyrs and an assumed reddening of E{B - V) = 0.05 
yield the best fit for our isochrones and are in good agree- 
ment with the latest edition of the Harris catalog ( Harris 1996[ l 
which gives [Fe/H] = -1.1 and - V) = 0.05. We further 
adopt [ff/Fe] = +0.4 fo llowing the observational measure- 
ments from Ful lton & Carney| ( jl996[ ), who found -hO.42 for 



the alpha elements in this cluster. 

The result is shown in Figure TBlfor the colours B-V, B-I, 
B - J, and B - K. With an adopted reddening of E{B - V) - 
0.05 and a standard reddening law with Ry - 3.1 we find an 
excellent match between the observed CMDs and models for 
all filter combinations. This test demonstrates that, at least 
in the intermediate metallicity regime, our isochrones fit data 
consistently well in all filter combinations in the case where 
there is no significant reddening. 

5.2.2. Isochrone Parameters for M4 

Incorrect assumptions for the isochrone parameters [Fe/H], 
[ff/Fe] or helium mass fraction (F) will inevitably lead to sys- 
tematic errors in the determination of Ry. For M4 we are us- 
ing isochi-ones with [Fe/H] = -1.0, [a/Fe] = +0.4, Y = 0.25 
and an age of 12 Gyrs which we will now justify. 

Although the main-sequence fiducial point which we are 
using to determine the colour excess is independent of age, 
it is important to know the metallicity of the cluster as accu- 
rately as possible since this parameter affects the shape and lo- 
cation of the isochrone, the absolute magnitude of the ZAHB, 



16 



Hendricks et al. 



TABLE 3 

Results for the individual filter combinations. 




[Fe/H]= -0.8 
[Fe/H]= -1.1 
[Fe/H]= -1.4 



3.5 4.0 
V-K 



5.0 



Fig. 17. — The overall shape of the the observed fiducial sequence is used to 
constrain the metallicity of M4 by isochrone fitting. Using 12 Gyr isochrones, 
we find a best fit when a value of [Fe/H] « -1.0 is adopted. 



and also the transformation-factor zero points from the red- 
dening law. There are many studies of the metallicity of M4, 
yielding values between [Fe/H] » -1.4 and -1.0 with recent 
results tending to favour the metal rich end of this interval. For 
example, jMucci arelli et al. (201 1 ) found [Fe/ H] ^ - 1.1 ±0.07 
for 87 MS and RGB stai's in M4 and Marino etal] ( |2008j de- 
rived [Fe/H] - -1.07 + 0.01 in their spectroscopic study of 
105 RGB stars. When we fit isochrones to our data for M4, 
we also find a best fit for a value of [Fe/H] close to -1.0. 
Figure [17] shows a comparison between isochrone fits with 
metallicities between [Fe/H] = -1.4 and -0.8, and shows that 
the overall slope of the observed fiducial sequence is matched 
best with [F e/H] - -1.0 and an a ge of 12 Gyrs. The re- 
cent work of |Bedin et iL] ( |2009| ) and |Mucciarelli et al.| ( f201 l| l 
supports this parameter choice: they both find a best fit for 
isochrones with [Fe/H] = -1.01 and age of 12 Gyrs derived 
from the WD cooling sequence. This latter nu mber is in good 
agreement with the results from Hansen et al.| ( |2004) who find 
a best fit for 12.1 Gyrs with the same method, which gives us 
confidence that our choice of isochrone parameters does not 
add a significant systematic error to the derived value of Ry- 
Since metallicity seems to be the parameter with the highest 
uncertainty and the biggest impact on Ry, we derive a quan- 
titative estimati on o f criRy) due to uncertain [Fe/H] assump- 
tions in section 1531 



Al - A2 
(1) 



Ofi'set 
(2) 



(3) 



Rv 

(4) 



For M4, r ecent alpha-element measurements from Marino 
eTaT] ([2008) yield [Mg/Fe] = h-0.50, [Ca/Fe] = -hO.28, 



TSiTFe] = H-0.48 and [Ti/Fe] = -hO.32 and therefore confirm 



B~V 0.37 ±0.01 1.00 

V-I 0.53 ±0.01 1.45 ±0.03 3.52 ±0.13 

V~J 0.97 ±0.02 2.63 ±0.06 3.64 ±0.11 

V-K 1.21 ±0.04 3.28 ±0.10 3.60 ±0.11 

B-I 0.90 ±0.02 2.45 ±0.07 3.52 ±0.28 

B'J 1.34 ±0.03 3.64 ±0.09 3.65 ±0.15 

B-K 1.57 ±0.05 4.28 ±0.12 3.60 ±0.14 

J-K 0.24 ± 0.03 0.65 ±0.07 3.52 ±0.28 

I- J 0.44 ±0.02 1.20 ±0.04 3.75 ±0.12 

I-K 0.68 ±0.03 1.84 ±0.08 3.66 ±0.13 

Note. — (2): Empirical colour offsets together with observational uncer- 
tainties between our data and a theoretical isochrone with [Fe/H] = -1.0, 
[t!'/Fe]= +0.4 and an age of 12 Gyrs. In (3) are the resulting transformation 
factors. (4) gives the individual best-fitting dust-type parameters for the line 
of sight to M4. The uncertainties in (3) and (4) represent the propagation of 
the error intervals from (2). 

the generally enhanced alpha-elements in old, metal-poor GC 
stars and justify our choice of [a/Fe] = 0.40. 

GC helium abundances are (still) assumed to be similar or 
equal to the primordial helium abundance arising from Big 
Bang nucleosynthesis. Th is would imply a he lium abundance 
near a value ofY- 0.25 jCyburt et aH2008| l, at least for the 
vast majority of clusters. There is no evidence in the literature 
that M4 is an exception to that rule. 

5.3. Results 

For the colour excess in B-V we find ^(B-y) = 0.37+0.01. 
Every colour offset we are using in direct comparison to 
E{B — V) yields a best-fitting value together with an accept- 
able uncertainty interval for Ry that produces the necessary 
observed excess for this filter combination. From those in- 
dividual results we derive the final value of Ry with the in- 
dependent sample to gain the smallest possible random and 
systematic errors. 

5.3.1. Individual Filter Combinations and the Independent Sample 

In a first step, we want to find out whether, in general, it 
is possible to reach a consistent fit for all colours within their 
photometric uncertainty limits. For that purpose we super- 
impose the individual results by shading the acceptable Ry 
intervals for each colour with the same intensity and in such a 
way that overlapping areas produce stronger shades. The top 
panel of Figure[18]shows such an overlay for the independent 
colours V-I, I- J, J-K and B-K. Using this colour selection, 
we find best fitting values in the range of Ry - 3.52 - 3.75 
with only a small interval where the model matches the data 
consistently for all colours at 3.62 < Ry < 3.65. 

The rather symmetric distribution around this best fitting 
region from colours where the filters alternate signs sug- 
gests that the deviations between different colours are mainly 
caused by systematic errors in individual filters. Similarly, it 
demonstrates that those errors are small compared to the ob- 
servational uncertainty, since a region of common overlap still 
exists. 

Next, we find the best fitting dust type for M4 by evaluat- 
ing the four selected colour combinations mentioned above, 
assuming that each of them yields an independent measure- 
ment of Ry. We determine the best fitting value of Ry and 
its uncertainty with the weighted mean and standard devia- 
tion of the results for the different filter combinations and find 
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Rv - 3.62 + 0.07. This result is in good agreement with pre- 
vious estimates in the Uterature, but with a total uncertainty 
several times smaller. 

The individual best fitting Ry values for the various colours 
as well as the measured colour offsets and the values of 
that are needed to produce these offsets are listed in Table [3] 
together with their respective uncertainties. Note that the 
value derived as the weighted mean nearly reproduces the re- 
gion of total overlap. Furthermore, all filter combinations de- 
duced from the independent sample yield best fitting values 
close to and symmetrically distributed around this mean. In 
Figure 19 we demonstrate the improvement in the isochrone 
fits when we use our derived value of Ry instead of a standard 
law with Rv = 3.1 and Table |4] we compare our result for Ry 
and its uncertainty to other results obtained in the past. Note 
that indirect measurements (i.e., along nearby lines of sight) 
can give a general hint concerning the dust type of M4, but 
cannot actually be used for an accurate prediction. 

It is important to point out that the dust- type parameter Ry 
determined here is not exactly the same as the actual total-to- 
selective extinction ratio AvUAb - Ay) for our filters. This is 
due to the fact that the dust type parameter Ry in the Cardelli 
equations is defined for Johnson's B and V bands using O 
and B type stars. As soon as a different filter set is used, Ry 
describes only the shape of the reddening law and not the ex- 
tinction ratio. The same is true if the reddening law is used for 
stars with significantly different temperatures, as is the case in 
our study. The difference is caused by the change in effective 
wavelengths for the B and V filters which is, in a sense, equiv- 
alent to using a different filter set. 

For the reddening correction and the distance determi- 
nation, however, we must use the actual ratio between 
Ay and E(B - V). The extinction values we obtain with 
a reddening law using Ay - 1.317 and Ry - 3.62 as 
applied to a star with the atmospheric parameters of our 
main-sequence fiducial point are listed in Table |5] and yield 
Av/{Ab - Ay) = 3.76 ± 0.07 where the confidence interval 
represents the uncertainty of the total-to-selective extinction 
ratio for a reddening law appropriate to these stars observed 
with this equipment. 
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Fig. 18. — Visualization of the determination of the dust type in the line of 
sight to M4: grey lines represent the dependence of F\^-\^ on Ry. Solid 
black horizontal lines indicate the best fitting value of F\^-x^ for each filter 
combination from observations and dotted black lines mark the uncertainty 
interval. Each filter combination produces a shaded region where the bound- 
aries are defined by the intersection of the dotted black lines with the cor- 
responding grey line. In the top panel, the results for the four independent 
combinations V-I,I-J,J-K and B - K are over-plotted in a way that 
overlapping areas produce a more saturated shade. 



TABLE 5 

Best fitting filter extinctions for M4 



Rv 


Au 
Av 


Av 


Ar 
Av 


A, 
Av 


Aj 
Av 


Ah 

Av 


Ak 
Av 


3.62 


1.414- 


1.266 


0.831 


0.608 


0.302 


0.191 


0.123 



TABLE 4 

Comparison of the determined value of Ry to ltterature values. 



Author 


Quoted Ry 


O" 


Method 


Notes 


Cudwoith & Rees (1991) 


3.3 


0.7 


direct 




Dixon & Longmore (1993) 


~ 4 




direct 




Peterson et al. (1995) 


~ 4 




direct 




Ivans et al. (1999) 


3.4 


0.4 


direct 




Chini(1981) 


4.2 


0.2 


indirect 


(3) 


Clayton & Cardelli (1988) 


3.8 




indirect 


(2), (3) 


Vrba(1993) 


3.99 


0.18 


indirect 


(1), (3) 


This Work 


3.62 


0.07 


direct 





Note. — (1): This result was obtained by the measurement of 12 stars in 
the Sc o-Oph cloud complex. The stated value (that is used by [Richer e t al.| 
(1997) and others) only represents the mean value of all measurements, where 
the individual results actually vary between Ry = 3.15andR\/ = 5.25(!). This 
study is an excellent example of the fact that it is not possible to infer the value 
of Ry for an object from its direct environment. Specifically, even if the lines 
of sight are similar, we cannot evaluate possible inhomogeneities transverse 
to the line of sight, or along the line of sight between the nearer and more 
distant objects. (2): This result is obtained from the measurement of only one 
star. (3) The value actually represents the observed ratio Av/{Ab-Av) rather 
than the dust-type parameter for the Cardelli et al. reddening law. 



5.3.2. A Quick Test 

The two colours V - J and B - I provide a special case 
since the values of their transformation factors cross near the 
canonical dust type at Ry » 3.2 with E{V - J) being smaller 
for Ry < 3.2 and E(B - I) being smaller for Ry > 3.2. Ad- 
ditionally, both filter combinations have about the same the- 
oretical colour around the MSTO for all metallicities below 
[Fe/H] a; -0.8. Therefore those colours provide a quick qual- 
itative test for the dust type: if both B-I and V -J have about 
the same colour around the MSTO in an observed CMD, a 
standard dust type can be assumed. If V - 7 is significantly 
redder than B - I, Ry > 3.2; if V - 7 is significantly bluer, 
Ry < 3.2. 

5.4. Sources of Systematic Uncertainties 

Since we are using observational data together with model 
assumptions, there are several different sources of uncertainty 
contributing to the total uncertainty of Ry. After we estimated 
the random error in section [5. 1.2| we now discuss the differ- 
ent possible sources of systematic uncertainty and their sig- 
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nificance. 

The assumption of [Fe/H] affects the position and shape of 
the isochrone, the magnitude of the ZAHB and also the zero 
points for the reddening law. Since the uncertainty in metal- 
licity appears to have the biggest impact on the accuracy of 
our derived Ry and consequently on the derived distance to 
M4, we estimate the change in Ry arising from an uncertainty 
in [Fe/H], by performing the whole procedure to derive Ry 
while varying [Fe/H] by ±0.2 dex. This includes the use of 
a different isochrone to determine the ZAHB luminosity and 
the colour offsets as well as a different set of zero points for 
the transformation factors in our reddening law. The result for 
[Fe/H] = -1.2 yields a best-fitting value for Ry = 3.56 ± 0.07 
with a predicted reddening of y) = 0.39. An assumption 
of [Fe/H] = -0.8 yields = 3.96+0.09 and = 0.30. 

By comparing these numbers to the result we got assuming 
[Fe/H] - -1.0, the systematic effect on Ry due to uncer- 
tainties in [Fe/H] is highly asymmetrical and shows much 
larger deviations towards higher metallicities. However, in 
our case, the high-metallicity scenario can clearly be ruled 
out not only from isochrone fitting and from the unreasonably 
low amount of reddening resulting from this assumption, but 
also from recent spectroscopic studies. Therefore, we allow 
only the lower metallicity as a possible scenario and conclude 
that an uncertainty of cr[Fe/H] - +0.2 dex yields a change of 
Ry ~ ±0.06. 

The value of Ry that best fits the data further depends upon 
the exact reddening law we are using, and explicitly on the 
zero points for the transformation factors Fa,-a-,- These val- 
ues differ between different literature sources and our own 
calculations, and the same would be the case for the value 
of Ry needed to match isochrones with observations. We can 
only estimate the uncertainty introduced by the definition of 
the reddening law zero points by assessing the change in Ry 
when we use different sets of values from the literature and 
our own. The effect on Ry is, in any case, smaller than 0.1 
which might define an upper limit for the uncertainty due to 
the definition of the reddening law. However, most of the dis- 
crepancy here is introduced by different assumptions for the 
stellar atmospheric conditions. Once those are assumed cor- 
rectly, the actual uncertainty arising from the zero points is 
only caused by the quality of the spectra and might be signif- 
icantly smaller than the number stated above. 

When we determine observational colour excesses in dif- 
ferent filter combinations relative to theoretical isochrones, 
we assume that the models actually predict the unreddened 
CMD colours correctly. If this is not the case, for example 
due to faulty CT-relations or Teff-scales, the r esult of Ry is bi- 
ased. We discussed this issue in section l572l where we found 
that the isochrones provide a consistent fit to the nearly unred- 
dened cluster NGC 6723. However, minor discrepancies for 
at least some individual filters are likely to exist as our results 
for the individual filter combinations show. 

The observed data are calibrated to the |Landolt| ( |1992| l 
standard system for Johnson-Cousins UBV(RI)c and to the 
2MASS standard system for the NIR counterparts J and K,. 
The accuracy of the match between our adopted instrumental 
bandpasses and those from the respective standard systems is 
determined by the level to which any systematic differences 
between both can be corrected. Inaccuracies in the calibra- 
tion to the standard system have an impact similar to that of 
faulty CT-relations as they systematically bias the colour off- 
set between model and data. Generally, the calibration for our 
data shows no systematic deviation from this standard sys- 
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Fig. 19. — Direct comparison of two different laws: In each panel, the dashed 
line shows an isochrone corrected with a standard reddening law (Ry = 3.1), 
while a solid line shows an isochrone if a reddening law with Ry = 3.62 is 
used instead. Since both connections are normalized to B-V, both isochrones 
have identical loci here. 

tem. However, a final intrinsic uncertainty o f as m uch as 
cr(colour) x ±0.03 (see discussion in Stetsonl|2005|l has to 
be considered, especially if the data were obtained in just one 
observing run, as is the case for our NIR photometry. The 
calibration of the optical photometry of M4 was derived in- 
dependently on 20 photometric nights distributed among 9 
observing runs, so the overall standard calibration should be 
significantly better than 0.01 mag in this case. 

By using alternating colours such as V - I, I - J, J - K, 
and B - K to determine the value of Ry, we largely eliminate 
any systematic error which is specific for the individual fil- 
ters / and J. For this reason, the uncertainties in CT-relations 
and the calibration to the standard system have only a minor 
impact on the quality of our result. However, the uncertainty 
in our knowledge of [Fe/H] and the exact reddening law zero 
point is not filter-specific and therefore has to be taken into 
account when the total uncertainty of Ry is estimated. 

5.5. Implications for the Distance Modulus ofM4 

The absolute distance to a globular cluster, or its true dis- 
tance modulus Z)M(), strongly depends on the reddening law 
and the assumed dust type, if it is derived from standard can- 
dles of a stellar population such as variable stars or the appar- 
ent magnitude of the ZAHB. In this case 

DMo = DMapp-AA (6) 
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and 



d[kpc] 



jO([£"Wo/5]+l) 

fooo 



(7) 



where DM^pp is the observed or apparent distance modulus 
and Aa the extinction correction for the filter A. However, the 
extinction is usually determined with the colour excess and 
expressed by E(B - V). Aa is then obtained with the ratio of 
absolute to relative extinction R^ '. 



Aa = RaE{B - V). 



(8) 



Each Aa is therefore unambiguously correlated with the dust 
type, as is Ry, and as a consequence DMq can only be deter- 
mined as a function of that parameter. 

We determine the apparent distance modulus by comparing 
the ZAHB V-band magnitude to models and find (m - M)v - 
12.66 ± 0.03, which is slightly smaller than most of the v alues 
we find in the literature, especially the value given in |Har-| 
^ 7996) 1 of (m 



ris 



12.82. This is partly because we 



are usmg the findings from the treatment of spatial differen- 
tial reddening in this cluster, which defines a slightly brighter 
zero-age horizontal branch than the initial, uncorrected data. 
The use of a slightly higher metallicity for the isochrones than 
assumed in Harris additionally decreases the distance modu- 
lus. 

Using a reddening law with Ry - 3.62 ± 0.09 (including 
cr[Fe/H]) and a reddening of E(B - V) = 0.37 + 0.01 de- 
termined with isochrone fitting, we get a total-to-selective- 
extinction of A(//(Ab - Ay) - 3.76 for our filters. The 
distance to M4 is finally calculated with Eq. (jTj), yielding 
c/ = 1.80 ± 0.05 kpc or (m - M)q = 1 1.28 ± 0.06. 

The stated uncertainty for the absolute distance of M4 in- 
cludes the observational uncertainty for Ry {cr{Ry) - 0.09; 
including a cr[Fe/H] of ±0.2), the observational uncertainty 
for the apparent distance modulus (o-(DMapj,) = 0.03), and 
the uncertainty for the reddening (o-(E{B - V) - 0.01). 

Some publications already use a higher val ue of Ry, and 
found a distance to M4 o f d = 1.72 + 0.14 kpc ( Hansen et al. 
2004 Richer et al.|1997| l. Compared to these estimates, our 
distance is slightly larger but still within their observational 
uncertainty (see Table |6]). 



TABLE 6 

Overview of recently derived absolute distances to M4. 



Author 



Distance 



o- Av/E(B-V) Notes 



This work 1.80 kpc 0.05 3.76 

Harris (2010 edition) 2.20 kpc - 3.20 (1) 

Bedin(2009) 1.86 kpc - 3.80 

Hansen (2004) 1.72 kpc 0.18 3.80 (2) 

Peterson (1995) 1.72 kpc 0.14 - (3) 

Note. — (1): Ry comes from |Cudworth & Rees|(T990) , which is the ref- 
erence for H arris distance estimate. (2): This result is identical to [Richer] 
|et al.|(T997) . The uncertainty does not include uncertainties in Ry. (3): De- 
rived irom astrometric measurements; independent of photometric standard 
candles. 



6. SUMMARY AND DISCUSSION 

In this work we have investigated the dust properties in the 
direction to the globular cluster M4 using near-infrared J and 
Ks observations obtained with SOFl on the New Technology 
Telescope in Chile, which have been combined with archival 
optical Johnson-Cousins UBV(RI)c photometry. We deter- 
mined the size and the distribution of the spatial differential 



reddening across the cluster face, the mean absolute amount 
of reddening, and finally the type of dust in the line of sight 
to M4. With this information we were able to determine the 
absolute distance to M4 with a higher precision than previous 
estimates. 

We have used colour excesses derived from a combination 
of NIR and optical photometry together with Victoria-Regina 
isochrones to determine the dust type parameter represented 
by Ry in the Cardelli reddening law (Cardelli et al. 1989) and 
the actual value of the total-to-selective extinction ratio for 
M4. Our method is independent of age assumptions, insensi- 
tive to metallicity and appears to be significantly more precise 
and accurate than existing spectroscopic approaches. In the 
following we summarize the main results: 

• We investigated the reddening variations across the face 
of the field of M4, finding a total peak-to-peak differ- 
ence of 6E(B - V) ~ 0.2 mag from the north-east to 
the south-west within a radius of 10' around the clus- 
ter centre. Therefore the differential effects within the 
field of M4 are about half the size of its mean total red- 
dening, which we found to be E(B - V) = 0.37 ± 0.01. 
By correcting for the variations in reddening, we have 
been able to decrease the observational scatter in our 
photometry by about 50% and consequently increase 
the precision of important parameters like the zero-age 
horizontal branch luminosity by the same factor 

• We have written a program to investigate the impact 
of stellar temperature, surface gravity and metallicity 
on the extinction properties in different broadband fil- 
ters and on the necessary reddening corrections. We 
found similar sized effects for temperature and surface 
gravity within typical globular cluster parameter lim- 
its; each causes a change of about 3% in the necessary 
correction factor for each filter combination. However, 
the impact of gravity seems to be significant only for 
the coolest main sequence and red giant branch stars. 
It is important to note that metallicity causes similar 
changes in the correction factors as the other atmo- 
spheric parameters when it is changed from solar metal- 
Ucity ([Fe/H] = 0.0) to [Fe/H] = -2.5. 

• We made use of three different spectral databases to in- 
tercompare the results obtained with synthetic ATLAS9 
spectra and observed spectra for different temperatures 
(using the atlas of Pickles ( 1998 )) and met allicities (us- 
ing the atlas of ^Sanchez-Blazquez et al.| ( 2006) ). We 
found excellent agreement between the difi'erent types 
of spectra for all metallicities and for temperatures 
above 5 000 K; however significant deviations are found 
at lower temperatures. Inadequate treatment of molecu- 
lar bands within the ATLAS9 models is suggested to be 
the most likely cause of the discrepancies at low tem- 
peratures. 

• We have used a combination of optical and near- 
infrared colours to determine the average dust type 
toward M4, as characterized by Ry in the Cardelli 
et al. reddening law, by fitting the newest generation of 
Victoria-Regina isochrones to our observed CMD, re- 
quiring a consistent fit for all colours. This method is 
independent of the assumed age of the system and in- 
sensitive to the assumed metallicity. With our method, 
we have determined the dust-type parameter to be Ry - 
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3.62 + 0.07 ± 0.06, where the first uncertainty is the 
observational precision and the second represents pos- 
sible systematic effects due primarily to uncertainties in 
[Fe/H] . Using this dust-type parameter, we find an ac- 
tual ratio for the total-to-selective extinction in our fil- 
ters of Av/{Ab - Ay) = 3.76. A generally large number 
for Rv has been proposed several times in the literature 
but never with such a strong quantitative justification. 

• We have tested our isochrones for optical and NIR 
colours in the intermediate metal poor regime on the 
nearly unreddened globular cluster NGC 6723 and con- 
firm that they provide a good and consistent fit for all 
filter combinations. It is worth noting that NGC 6723 
is identical in age and metallicity to M4 within the pho- 
tometric uncertainties, with a similar bimodal HB se- 
quence. It is therefore an essentially unreddened twin 
to the highly reddened M4. 

• With the new dust type and the knowledge about the 
spatial differential reddening effects, we have reexam- 
ined the distance to M4 and find a significantly (~ 18%) 
smaller value for the a bsolute distanc e compared to the 
parameter provided by Harris ( 1996 1. This makes M4 
unambiguously the closest globular cluster to the Sun 
with a distance of 1.80 ± 0.05 kpc. Our value is, how- 
ever, somewhat larger than some recent estimates of 
d = 1.72kpc, based on TJy = 3.8 ( .Hansen et al.,2004j ). 

The proper treatment of interstellar extinction is a crucial 
step to achieve accurate astrophysical results from photomet- 
ric observations. As seen from this perspective, M4 is an ex- 
cellent object to study the effect of interstellar extinction and 
its impact on the accuracy of the derived parameters due to the 
strong, differential and unusual extinction it suffers in com- 
bination with its extreme proximity, which allows extremely 
deep and precise photometry in all photometric filters. In a 
way, it serves as an in-situ laboratory from which we can learn 
about the application of extinction corrections for fainter ob- 



jects such as stellar populations in the bulge of our Galaxy or 
resolved extragalactic populations. 

Our study reveals that individual — object-specific — 
parameters of the target of investigations such as temperature, 
surface gravity and metallicity can all have a significant 
impact on the extinction properties of broadband filters 
and consequently on the colour and extinction correction 
procedure. Even though variations of these parameters within 
the population do not seem to add a significant systematic 
error for systems with E{B - V) < 0.5, the systematic 
difference caused by the use of wrong zero points (e.g., 
defined for a Vega-like star) can be much larger and should 
be taken into account for all reddened systems. Therefore 
we recommend the use of an object-specific reddening law 
defined for appropriate values of the crucial target parameters 
for the correction procedure. However, For highly reddened 
targets with E(B - V) = 0.5 or more, the use of a star-by-star 
correction should be considered to account for the differential 
effects within the population. 

Moreover, our work on M4 has shown that the shape of 
the reddening law defined by the type of dust in the line of 
sight can be as object-specific as the stellar parameters. Espe- 
cially for highly reddened objects, the unconsidered use of the 
standard dust-type parameter for the interstellar medium — or 
other generalized reddening laws — can cause significant sys- 
tematic errors in the interpretation of photometric data. 

Finally, we find that model atmospheres and the associated 
synthetic spectra still have problems adequately predicting the 
complex properties of molecular bands at cool temperatures. 
In terms of the derived reddening-correction factors, this ef- 
fect is small. However, more work needs to be done to investi- 
gate the impact of this insufficiency on the actual CT-relations 
that are used to compare model predictions to observations. 

B. Hendricks gives thanks to Karsten Brogaard for many 
helpful discussions of several reddening-related issues which 
are incorporated in this manuscript. This work was supported 
by the Natural Sciences and Engineering Research Council of 
Canada through a Discovery Grant to DAV. 
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APPENDIX 



FIDUCIAL SEQUENCES 



TABLE 7 
Empirical fiducial sequence for M4. 



V 


(y-7) 


iV-K) 


Vc 


iV-J)c 


(y-K)c 


12.550 


2.952 


3.922 


12.408 


2.930 


3.884 


12.850 


2.840 


3.758 


12.715 


2.819 


3.721 


13.150 


2.771 


3.670 


13.014 


2.750 


3.633 


13.450 


2.718 


3.605 


13.314 


2.697 


3.568 


13.750 


2.658 


3.533 


13.620 


2.638 


3.498 


14.050 


2.618 


3.463 


13.925 


2.601 


3.429 


14.350 


2.572 


3.397 


14.246 


2.560 


3.369 


14.650 


2.537 


3.336 


14.566 


2.527 


3.313 


14.950 


2.510 


3.279 


14.870 


2.500 


3.258 


15.250 


2.478 


3.244 


15.176 


2.471 


3.224 


15.550 


2.459 


3.213 


15.488 


2.453 


3.197 


15.850 


2.425 


3.170 


15.777 


2.423 


3.150 


16.150 


2.397 


3.126 


16.085 


2.398 


3.108 


16.300 


2.344 


3.046 


16.249 


2.343 


3.032 


16.365 


2.260 


2.928 


16.326 


2.261 


2.918 


16.450 


2.155 


2.758 


16.427 


2.157 


2.752 


16.600 


2.108 


2.675 


16.597 


2.112 


2.674 


16.750 


2.082 


2.650 


16.759 


2.087 


2.652 


17.050 


2.076 


2.647 


17.066 


2.089 


2.651 


17.350 


2.094 


2.674 


17.359 


2.106 


2.676 


17.650 


2.130 


2.727 


17.650 


2.140 


2.727 


17.950 


2.176 


2.802 


17.949 


2.185 


2.802 


18.250 


2.245 


2.902 


18.247 


2.253 


2.901 


18.550 


2.324 


3.021 


18.546 


2.331 


3.020 


18.850 


2.405 


3.138 


18.838 


2.408 


3.135 


19.150 


2.516 


3.308 


19.103 


2.513 


3.296 


19.450 


2.627 


3.460 


19.383 


2.620 


3.442 


19.750 


2.733 


3.619 


19.707 


2.731 


3.608 


20.050 


2.870 


3.806 


19.965 


2.863 


3.783 


20.350 


2.995 


3.978 


20.264 


2.985 


3.954 


20.650 


3.120 


4.140 


20.568 


3.112 


4.118 


20.950 


3.248 


4.297 


20.862 


3.240 


4.274 


21.250 


3.341 


4.408 


21.161 


3.333 


4.384 


21.550 


3.446 


4.505 


21.459 


3.437 


4.480 


21.850 


3.564 


4.623 


21.758 


3.554 


4.598 


22.150 


3.651 


4.704 


22.056 


3.640 


4.679 



NoFE. — V- J and V - are the fiducial lines derived without any correction. For (V - J)c and (V - K)^ a star-by-star correction has been applied appropriate 
to the atmospheric and chemical parameters of the individual stars and normalized to MSTO conditions {T^g = 6 000 K, log g = 4.5). 
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TABLES 

Empirical fiducial sequence for NGC6723. 



V 


(V - J) 


{V - K) 


12.850 


2.878 


3.861 


13.150 


2.635 


3.593 


13.450 


2.473 


3.368 


13.750 


2.348 


3.170 


14.050 


2.271 


3.042 


14.350 


2.176 


2.904 


14.650 


2.061 


2.771 


14.950 


2.017 


2.684 


15.250 


1.959 


2.596 


15.550 


1.887 


2.503 


15.850 


1.836 


2.429 


16.150 


1.779 


2.365 


16.450 


1.757 


2.317 


16.750 


1.716 


2.264 


17.050 


1.691 


2.221 


17.350 


1.673 


2.195 


17.650 


1.636 


2.155 


17.950 


1.620 


2.112 


18.250 


1.567 


2.037 


18.400 


1.439 


1.866 


18.550 


1.335 


1.679 


18.850 


1.255 


1.592 


19.150 


1.249 


1.577 


19.450 


1.267 


1.591 


19.750 


1.305 


1.652 


20.050 


1.331 


1.714 


20.350 


1.419 


1.824 


20.650 


1.446 


1.931 


20.950 


1.529 


2.050 


21.250 


1.665 


2.248 


21.550 


1.831 


2.394 


21.850 


2.058 


2.695 


22.150 


2.220 


2.867 


22.450 


2.381 


3.053 



TABLE 9 

Empirical ZAHB fiducial sequence for M4. 



(V-J) jV-K) V 



1.15 


1.35 


13.78 


1.23 


1.46 


13.58 


1.28 


1.55 


13.52 


1.34 


1.71 


13.49 


1.42 


1.85 


13.50 


2.04 


2.60 


13.46 


2.24 


2.84 


13.47 


2.29 


2.92 


13.47 


2.40 


3.09 


13.40 
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TABLE 10 

Empirical ZAHB fiducial sequence for NGC 6723. 



{V -J) {V-K) V 



-0.04 


-0.05 


16.81 


0.05 


0.05 


16.42 


0.18 


0.20 


15.88 


0.31 


0.30 


15.61 


0.41 


0.50 


15.53 


0.60 


0.76 


15.50 


1.13 


1.46 


15.53 


1.49 


1.87 


15.53 


1.62 


2.13 


15.33 



DIFFERENTIAL EXTINCTION TABLES 



TABLE 1 1 

Extinction coefficients for different temperatures. 



log Teff 


Au 
Av 


As 
Av 


Ar 
Av 


A; 
Av 


Aj 
Av 


Ah 
Av 


Ak 
Av 


Aviabs) 


4.600 


1.578 


1.331 


0.832 


0.580 


0.284 


0.178 


0.115 


1.147 


4.450 


1.581 


1.330 


0.831 


0.580 


0.283 


0.178 


0.115 


1.146 


4.280 


1.572 


1.325 


0.828 


0.579 


0.284 


0.178 


0.115 


1.145 


4.150 


1.564 


1.323 


0.828 


0.579 


0.284 


0.177 


0.115 


1.144 


4.070 


1.559 


1.324 


0.829 


0.579 


0.286 


0.177 


0.116 


1.143 


4.030 


1.555 


1.321 


0.829 


0.579 


0.284 


0.178 


0.115 


1.142 


3.980 


1.552 


1.317 


0.828 


0.580 


0.284 


0.178 


0.115 


1.142 


3.929 


1.550 


1.317 


0.826 


0.579 


0.284 


0.178 


0.115 


1.140 


3.858 


1.557 


1.315 


0.823 


0.581 


0.285 


0.179 


0.116 


1.137 


3.815 


1.561 


1.314 


0.822 


0.580 


0.284 


0.179 


0.116 


1.135 


3.778 


1.560 


1.311 


0.820 


0.580 


0.285 


0.180 


0.116 


1.134 


3.764 


1.559 


1.310 


0.819 


0.580 


0.285 


0.180 


0.116 


1.133 


3.747 


1.560 


1.306 


0.818 


0.581 


0.285 


0.180 


0.116 


1.132 


3.715 


1.559 


1.305 


0.817 


0.582 


0.285 


0.180 


0.117 


1.129 


3.689 


1.567 


1.300 


0.818 


0.582 


0.285 


0.180 


0.117 


1.125 


3.653 


1.565 


1.303 


0.817 


0.584 


0.286 


0.181 


0.117 


1.122 


3.638 


1.567 


1.298 


0.815 


0.584 


0.287 


0.181 


0.118 


1.120 


3.622 


1.569 


1.299 


0.817 


0.586 


0.288 


0.182 


0.118 


1.116 


3.602 


1.568 


1.299 


0.810 


0.585 


0.289 


0.182 


0.118 


1.114 


3.580 


1.565 


1.298 


0.805 


0.583 


0.287 


0.182 


0.118 


1.115 


3.566 


1.557 


1.296 


0.797 


0.582 


0.288 


0.181 


0.118 


1.116 


3.550 


1.556 


1.299 


0.792 


0.581 


0.287 


0.182 


0.117 


1.116 


3.519 


1.548 


1.295 


0.771 


0.576 


0.286 


0.181 


0.117 


1.119 


3.470 


1.550 


1.285 


0.739 


0.571 


0.284 


0.181 


0.116 


1.120 




1.664 


1.321 


0.819 


0.594 


0.276 


0.176 


0.112 


S98 




1.569 


1.323 


0.751 


0.478 


0.282 


0.192 


0.116 


C89 






1.32 


0.82 


0.59 


0.29 


0.23 


0.11 


B98 




1.545 


1.326 


0.795 


0.558 


0.267 


0.169 


0.114 


M04 



Note. — Relative extinctions for Johnson-Cousins UBV(RI)c (Landolt 1992i and 2MASS J, H, and filters for spectra with Ry = 3.10, logg = 4.5, 
[Fe/H] = 0.0 and E(B - V) = 0.36 for dift'erent temperatures. Since the ArLAS9 models c annot make accurate predicti ons for temperatures below 5 000 K, 
observ ed fluxes are used to calculate values in this table. S 98: [Schlegel et al.H1998f ; C 89: [CardeUi et aL](T989) ; B 98: [Bessell et al"1(T998t ; M 04: |McCall| 
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TABLE 12 

Extinction coefficients for different metallicities. 



[Fe/H] 


Av 

Ay 


An 
Av 


Ar 
Av 


A, 

Av 


Aj 
Av 


Ah 
Av 


Ak 
Av 


Av(abs) 




+0.5 


1.565 


1.308 


0.821 


0.582 


0.284 


0.180 


0.116 


1.132 




+0.2 


1.565 


1.310 


0.821 


0.582 


0.284 


0.180 


0.116 


1.132 




+0.0 


1.565 


1.311 


0.821 


0.581 


0.284 


0.180 


0.116 


1.133 




-0.5 


1.566 


1.313 


0.820 


0.581 


0.284 


0.180 


0.116 


1.133 




-1.0 


1.567 


1.316 


0.820 


0.580 


0.284 


0.180 


0.116 


1.134 




-1.5 


1.568 


1.318 


0.819 


0.580 


0.284 


0.179 


0.116 


1.134 




-2.0 


1.569 


1.319 


0.819 


0.580 


0.284 


0.179 


0.116 


1.134 




-2.5 


1.570 


1.320 


0.819 


0.580 


0.284 


0.179 


0.116 


1.134 




-4.0 


1.571 


1.321 


0.819 


0.580 


0.284 


0.179 


0.116 


1.134 






1.664 


1.321 


0.819 


0.594 


0.276 


0.176 


0.112 




S98 




1.569 


1.323 


0.751 


0.478 


0.282 


0.192 


0.116 




C89 






1.32 


0.82 


0.59 


0.29 


0.23 


0.11 




B98 




1.545 


1.326 


0.795 


0.558 


0.267 


0.169 


0.114 




M04 



Note. — Relative extinctions for Johnson-Cousins UBV(RI)c (Landolt 19921 and 2MASS J, H, and fi lters for spec tra with Teg = 6000K, log g = 4.5, 
Ry = 3.10 and E(B - V) = 0.36 for difl'erent metallicities. S 98: jSchlegel et al.,^1998^ ; C 89: [CardelU et al.|(l989( ; B98: |Bessell et al.|(l998/ ; M 04: |McCall| 



TABLE 13 

Extinction coefficients for different dust type parameter Ry. 



Rv 


Au 
Av 


Ab 
Av 


Ar 
Av 


A, 

Av 


Aj 
Av 


Ah 
Av 


Ak 
Av 


Av(ahs) 




2.6 


1.687 


1.369 


0.803 


0.545 


0.260 


0.165 


0.106 


0.953 




2.8 


1.633 


1.343 


0.811 


0.561 


0.271 


0.171 


0.110 


1.025 




3.0 


1.586 


1.321 


0.818 


0.575 


0.280 


0.177 


0.114 


1.097 




3.2 


1.545 


1.301 


0.824 


0.587 


0.288 


0.182 


0.118 


1.168 




3.4 


1.509 


1.284 


0.829 


0.598 


0.295 


0.187 


0.120 


1.240 




3.6 


1.477 


1.268 


0.834 


0.608 


0.301 


0.191 


0.123 


1.312 




3.8 


1.448 


1.255 


0.838 


0.616 


0.307 


0.194 


0.125 


1.384 




4.0 


1.421 


1.242 


0.842 


0.624 


0.312 


0.197 


0.127 


1.456 




4.2 


1.398 


1.231 


0.845 


0.631 


0.317 


0.200 


0.129 


1.527 






1.664 


1.321 


0.819 


0.594 


0.276 


0.176 


0.112 




S98 




1.569 


1.323 


0.751 


0.478 


0.282 


0.192 


0.116 




C89 






1.32 


0.82 


0.59 


0.29 


0.23 


0.11 




B98 




1.545 


1.326 


0.795 


0.558 


0.267 


0.169 


0.114 




M04 



Note. — Relative extinctions for Johnson-Cousins U BV(RI)c (La ndolt' 1992 1 and 2MA SS J, H and filter s for spectra w ith T^s = 6000K, log^ = 4.5, 
[Fe/H] = 0.0 and E(B - V) = 0.36 for difterent Ry. S 98: |Schlegel et al. 1 1998 1; C 89: Cardelli et al.|(T989) ; B 98: |Bessell et"arifl998 l; M 04: McCall. l 2004. 1. 



